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SUMMARY 


.  \ 

'  It  has  been  demonstrated  that  a  strong  interfaciai  bond  is  required  to 
produce  the  observed  strengthening  in  discontinuous  SiC/Al  composites. 
Dislocation  generation  (which  is  the  strengthening-mechanism)  is 
dependent  upon  a  good  bond  existing  between  the^SiCjand  the  (Ai /matrix. 

Although  it  has  been  shown  that  the  stresses  developed  in  the  matrix 
as  a  result  of  the  differences  in  the  coefficients  of  thermal  expansion 
between  SiC  and  Ai  are  greatly  reduced  by  the  motion  and  generation  of 
dislocations,  an  elastic  residual  stress  can  still  remain. 

An  investigation  was  undertaken  and  it  was  determined  that  in  the 
discontinuous  SiC/Al  composites  a  tensile-elastic-residual  stress 
remains  in  the  matrix.  This  determination  is  based  on  results  from  a 
theoretical  model,  tensile  and  compression  testing  and  X-ray  measure- 
m.ents.  The  results  from  all  the  tests  and  the  model  are  in  agreement. 
The  magnitude  of  the  residual  stress  is  small,  e.g.,  with  20  V%  whisker 
in  an  annealed  6061  AI  alloy  matrix,  the  residual  stress  is  70  MPa.  If 
the  matrix  is  annealed  1100  AI,  then  residual  stress  is  still  smaller: 
43  MPa. 

The  tensile  residual  stress  produces  a  Bauschinger  Effect  (BE)  which 
is  highly  dependent  upon  the  initial  direction  of  loading.  If  the  initial 
direction  of  loading  is  in  tension,  the  BE  is  small,  if  the  initial  direc¬ 
tion  of  loading  is  in  compression,  the  BE  is  very  large. 

It  has  been  observed  that  as  the  strength  of  the  matrix  increases, 
e.g.,  due  to  alloying  or  age  hardening,  the  percentage  increase  in  yield 


^or  ultimate  tensile  s^ength  decreases  for  a  given  volume  per  cent 
SiC/Al  composite.<^If  the  matrix  is  1100  Al  and  20  V%  SiC  whisker  is 
added,  then  there  is  a  420%  increase  in  strength,  whereas,  if  the  matri 
is  annealed  6061  Al  alloy,  there  is  a  ^  181%  increase  in  strength,  .^or 
annealed  7091  matrix  the  increase  is  .v  190%.  If  the  absolute  magnituc 
is  considered,  the  increase  is  greatest  for  1100  Al  matrix  and 
decreases  with  increasing  alloy  content.  The  basic  reason  for  this 
di.*^ference  in  strengthening  is  due  to  the  change  in  the  magnitude  of  the 
dislocation  density  and  of  the  tensile  residual  stress.  In  the  case  of  the 
1100  matrix,  there  is  maximum  dislocation  generation  due  to  di.^.^r- 
ences  in  coefficients  of  thermal  expansion  and  a  minimum  tensile 
residual  stress,  because  the  yield  stress  of  annealed  1100  matrix  is 
less  than  the  yield  stress  of  the  annealed  6061  and  7091  matrix. 

In  the  .‘‘irst  section  is  a  list  of  presentations,  and  the  second  section 
includes  a  list  of  papers  in  preparation.  The  third  section  is  a  com.pila 
tion  of  the  oubiications  which  have  been  oroduced  under  this  contract. 
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Introduction 


A  basic  premise  in  most  theories  of  composite  strengthening 
is  that  a  bond  of  some  type  exists  between  the  matrix  and  the 
reinforcement.  This  bond  is  Involved  in  the  mechanism  of  load 
transfer  between  the  matrix  and  the  reinforcement.  Therefore,  in 
terms  of  tensile  strength  a  good  or  sufficient  bond  is  required 
to  obtain  the  maximum  strength  from  the  composite.  In  the  model 
metal-matrix  composite,  i.e.,  the  Cu-W  composite,  the  bond 
strength  between  W  and  Cu  is  very  good^^^.  In  the  case  of  SiC-Al 
composites  there  have  been  no  direct  measurements  of  the  bond 
strength,  but  indirect  data  Indicates  that  the  bond  strength  is 
quite  good^^*^^. 

It  is  presumed  that  the  bond  strength  is  related  to  the 
nature  of  the  Interface  between  the  matrix  and  the  reinforcement. 
There  have  been  three  mechanisms  or  hypothesis  put  forward 
concerning  the  Interface  between  the  SiC  and  the  A1  matrix.  The 
first  is  that  there  is  a  SIO2  film  or  layer  between  the  SIC  and 
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the  A1  matrix,  this  is  an  idea  put  forward  by  Parratt^^^.  The 
second  hypothesis  is  that  there  is  nothing  at  the  interface. 

This  hypothesis  is  based  on  the  work  that  has  been  done  on  the 
model  metal-matrix  composite,  i.e.,  the  Cu-W  composite^^^ .  In 
that  particular  system  no  Interphases  were  found  between  the  Cu 
and  W.  The  third  hypothesis  is  that  there  is  a  Al/jC^  film 
between  the  SIC  and  the  A1  matrix,  this  concept  was  put  forward 
on  a  thermodynamics  basis.  If  molten  A1  exists  with  a  small 
amount  of  dissolved  Si,  then  it  thermodynamically  possible  for 
AlijC^  to  form^^\ 

The  purpose  of  this  Investigation  was  to  determine  in 
greater  detail  the  nature  of  the  interface  between  SiC  and  A1  and 
W  and  Cu. 

Results 

A  scanning  Auger  microprobe  analysis*  was  undertaken  of 
fractured  SiC-Al  samples,  in  particular  fractured  whisker  samples 
produced  by  ARCO  SILAG.  These  samples  were  fractured  in  the 
scanning  Auger  ralcroprobe  and  examined  shortly  thereafter.  When 
the  fracture  surfaces  were  examined  there  were  no  differences,  in 
signal  as  a  function  of  position,  in  other  words,  a  100%  A1 
signal  was  obtained  from  the  matrix  as  well  as  from  all  possible 
areas  on  the  SIC  whiskers  as  shown  in  Pig.  1.  A  possible 

•  The  Auger  microprobe  analysis  were  performed  at  the  Center  for 
Research  in  Surface  Science  and  Submicron  Analysis,  Montana 
State  University. 
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explanation  of  the  essentially  lOOj  A1  content  on  the  whisker  Is 
that  the  A1  could  just  be  mechanically  adhering,  very  similar  to 
the  butter  knife  anology.  (If  you  stick  a  butter  knife  Into 
butter  and  pull  It  out  you  have  a  coating  of  butter  on  the  butter 
knife.)  Therefore,  a  series  of  etching  experiments  were 
conducted  In  the  microprobe  and  It  was  only  after  considerable 
sputtering,  (which  Is  basically  an  Ar  Ion  milling  operation)  that 
the  SI  and  C  signals  were  obtained.  The  sputter  etching  was 
Interrupted  at  Intervals  and  the  surfaces  of  the  whiskers 
examined,  there  were  no  Indications  of  SIO2,  Al20^,  Al^O^*  This 
result  led  us  to  the  belief  that  the  A1  had  diffused  Into  the 
SIC.  In  order  to  determine,  whether  this  Is  a  real  effect,  and 
not  a  sputtering  artifact,  several  samples,  which  were  metalo- 
graphlcally  polished,  were  examined  In  the  Auger  microprobe.  The 

data  from  a  platelet  of  SIC  sample  from  DWA,  are  shown  In  Pig.  2. 

•• 

In  the  SEM  micrograph  the  scanning  Auger  probe  traces  are 
superimposed.  One  trace  Is  for  A1  and  other  trace  Is  for  SI. 

This  shows  that  A1  has  penetrated  or  diffused  Into  the  SIC,  but 
there  Is  no  SI  or  C  In  the  matrix.  Then  this  sample  was  sputter 
etched  and  the  process  repeated  and  basically  the  same  result  was 
obtained.  In  order  to  determine  what  the  effect  of  sputtering, 
(l.e.,  could  the  sample  surface  become  cross  contaminated  by  the 
sputtering  operation)  the  sample  was  removed,  etched  slightly 
with  Keller's  etch,  and  then  put  back  Into  the  microprobe  and 
then  the  process  was  repeated  and  again  the  same  results  were 
obtained.  Figure  3  Is  the  detailed  Auger  probe  analysis  of  the 
SI  and  A1  traces,  the  C  trace  superimposes  right  on  the  SI 
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trace.  The  interface  is  then  defined  as  having  finite  width  and 
that  is  the  maximum  probe  size  that  is  associated  with  the 
scanning  microprobe  apparatus. 

In  order  to  further  validate  these  results,  a  STEM  analysis 
was  conducted  of  SIC  samples,  and  Pig.  4  is  a  result  of  an  energy 
dispersive  X-ray  analysis  from  the  STEM  using  a  thin  foil  of  the 
composite.  It  is  obvious  that  again  tnere  is  evidence  of  A1 
diffusion  into  the  SiC.  In  this  particular  case  the  probe  size 
is  much  smaller  (~  50  A)  than  that  in  the  scanning  Auger 
analysis.  The  thickness  of  the  specimen  was  1000  -  1500  A; 
(accelerating  voltage,  of  the  electrons  200  KV)  thus  the  spatial 
resolution  of  the  technique  was  about  300  A,  several  times 
smaller  than  the  range  of  penetration  observed.  In  order  to 
determine  whether  this  may  be  a  universal  affect,  an  analysis  was 
begun  of  Cu-W.  Figure  5  is  a  typical  scanning  Auger  microprobe 
result  of  Cu  and  W.  The  data  Indicates  that  Cu  has  diffused  into 
the  W,  and  there  is  a  very  sharp  edge  of  the  W  filament  and  the 
matrix.  It  should  be  kept  in  mind  that  the  experimental  results 
of  the  bond  strength  between  Cu  and  W  is  extremely  good.  STEM 
analysis  of  Cu-W  interfaces  are  in  progress. 

Discussion  of  Conclusion 

The  observations  of  A1  penetration  into  SIC  and  the  Cu 
penetration  into  W  to  the  extent  observed  is  unexpected.  There 
have  been  at  least  two  Investigations  of  A1  diffusion  into 
Sic(7,8),  from  those  investigations  it  is  possible  to 
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calculate  the  mean  diffusion  distance  X  at  773  K, 


X  =  /Dt 


where  D  =  D^e  -  Q/RT,  t  is  the  time  of  diffusion,  Q  for  Al-Sic 

system  is  471  kJ/mole  and  equals  1.8  cm^/sec.  The  value  of  X 

is  less  than  one  A  per  year.  (The  values  of  and  Q  given  in 
literature  vary  substantially  from  each  other.  However  all  of 
these  values  give  a  value  of  X  orders  of  magnitude  lower  than  the 
observed  penetration.)  Needless  to  state  this  does  not  agree 

with  the  present  data.  In  the  case  of  Cu  diffusion  into  W  there 

have  been  no  known  reported  investigations.  In  fact  there  is  no 
phase  diagram  of  Cu-W^^\  but  there  is  SiC-Al^Cii  phase 
dlagrara^^®^  at  2133  K  and  it  shows  solubility  of  Al^C^i  in  SiC. 

The  fact  that  there  is  solubility  of  Al^C^  in  SiC  does  not 
explain  how  the  A1  could  diffuse  that  far.  However,  from  the 
second  law  of  thermodynamics  the  presence  of  A1  is  SiC  and  Cu  in 
W  can  be  explained,  as  due  a  reduction  in  free  energy  due  to  the 
entropy  of  mixing.  In  fact  there  is  some  data  (at  2073  K)  of  the 
solubility  of  Al^^^^  in  SIC.  The  depth  of  penetration,  l.e., 
diffusion,  cannot  be  explained  in  terms  of  classical  bulk 
diffusion,  there  must  be  a  short  circuit  path.  In  the  case  of 
Cu-W,  there  may  be  substantial  grain  boundary  diffusion,  due  to 
the  small  grain  size  in  the  W  filaments. 

In  terms  of  bonding  between  the  reinforcement  and  the 
matrix,  in  both  cases  the  bonding  appears  to  be  very  good,  and 
there  is  penetration  of  the  matrix  into  the  reinforcement.  It 


may  be  concluded  that  the  diffusion  bonding  results  in  a  high 
strength  bond. 
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Pl«.  1.  A  SEM  micrograph  taken  In  the  Auger  ralcroprobe  of  the 
fracture  surface  of  a  whisker  SlC-Al  composite. 


Pig.  2.  A  SEM  micrograph  taken  In  the  Auger  microprobe,  on  which 
the  A1  and  Si  Auger  probe  traces  have  been  superimposed. 
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A  STEM  EDXA  analysis  of  A1  concentration  as  a  function 
of  location  (I/^i/Isi  Plotted  In  arbitrary  units). 
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ABSTRACT 

Plastic  strains  and  the  extent  of  the  plastic 
zone  due  to  differential  thermal  expansion 
were  experimentally  determined  in  an  Al-SiC 
composite.  The  combined  plastic  shear 
strains  y^pss  the  Al-SiC  interface  for 
furnace-cooled,  air-cooled  and  quenched 
samples  were  found  to  be  1.32%,  1.23%  and 
0.99%  respectively.  Profiles  of  ^cpss  were 
plotted  versus  distance  from  the  interface  and 
compared  with  the  theoretical  distribution  of 
effective  strain  e.  The  theoretical  extent  of 
the  plastic  zone  measured  from  the  interface 
was  found  to  be  1.3  multiplied  by  the  particle 
radius.  This  value  was  slightly  less  than  the 
observed  value.  The  plastic  deformation  on 
the  heating  half  of  the  thermocycle  was  found 
to  be  at  least  equal  to  the  deformation  on  the 
cooling  half.  A  theoretical  treatment  of  the 
local  plastic  deformation  in  a  short  composite 
cylinder  was  suggested,  from  which  the  effec¬ 
tive  plastic  strain  e  and  the  extent  of  the 
plastic  zone  were  determined. 

1.  INTRODUCTION 

When  a  composite  material  is  subjected  to  a 
temperature  change,  local  plastic  deformation 
can  occur.  The  plastic  deformation  is  due  to  a 
stress  created  by  the  difference  between  the 
coefficients  of  thermal  expansion  of  the  com¬ 
ponent  phases,  and  this  stress  is  at  the  matrix- 
reinforcement  interfaces.  The  magnitude  of 
the  stress  is  equal  to  Aa  AT  where  Aa  is  the 
difference  between  the  coefficients  of  thermed 
expansion  of  the  phases  under  consideration 
and  AT  is  the  temperature  change.  The  rela¬ 
tively  large  (10  to  1)  difference  between  the 
coefficient  of  thermal  expansion  of  SiC  and 
that  of  aluminum  should  result  in  the  creation 
of  a  substantial  misfit  strain  at  the  Al-SiC 
interface  during  cooling  from  the  fabricating  or 
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annealing  temperature.  The  mechanical 
properties  of  the  composite  should  be 
affected  by  the  magnitude  and  extent  of  the 
plastic  deformation  that  takes  place  in  the  soft 
matrix  around  a  hard  particle  as  a  result  of 
the  misfit  relaxation  in  the  interface  region. 
Recent  investigations  [1-4]  have  indeed 
shown  the  important  role  of  the  interfacial 
regions  in  the  composite  strengthening.  Thus, 
knowledge  of  the  magnitude  of  plastic  strains 
and  the  size  of  the  plastic  zone  as  well  as 
other  characteristics  of  the  interfacial  region 
(the  bond  strength,  the  microstructure  etc.) 
should  contribute  to  the  understanding  of  the 
mechanism  of  composite  strengthening. 

The  experimental  determinations  of  plastic 
strains  and  plastic  zone  radii  about  a  particle 
in  the  matrix  due  to  the  difference  between 
the  coefficients  of  thermal  expansion  of  the 
phases,  to  our  knowledge,  have  not  been 
reported  in  the  literature. 

Several  theoretical  investigations  have  been 
undertaken  to  predict  the  magnitude  of  the 
plastic  strain  in  the  plastic  zone  around  a 
particle.  The  relaxation  of  the  misfit  caused 
by  the  introduction  of  an  oversized  spherical 
particle  into  a  spherical  hole  in  the  matrix 
has  been  analytically  described  by  Lee  et  al. 
[5].  Using  the  misfitting  sphere  model,  they 
calculated  strains  in  the  plastic  zone  that 
surrounds  a  hard  sphere  and  also  the  plastic 
zone  radius.  Hoffman  [6]  calculated  the 
overall  total  strains  in  the  tungsten -fiber- 
reinforced  80Ni-20Cr  matrix,  using  a  thick- 
wall  long-cylinder  approach  and  assuming 
that  a  hydrostatic  stress  state  exists  within 
each  constituent. 

Garmong  [7],  assuming  uniformity  of  the 
stresses  and  strains  in  the  matrix,  calculated 
deformation  parameters  for  a  hypothetical 
eutectic  composite  and  reported  values  of 
matrix  plastic  strains  that  were  of  the  order 
of  0.4%. 
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Dvorak  et  al.  [8]  developed  a  new  axi- 
symmetric  plasticity  theory  of  fiberous  com¬ 
posites  involving  large  thermal  changes.  The 
long-composite-cylinder  model  was  adopted 
as  a  composite  unit  cell,  and  the  microstress 
distribution  as  well  as  the  yielding  surfaces 
were  obtained  for  Al-W  composites. 

Mehan  [9]  calculated  the  residual  strains  in 
an  Al-a-Al203  composite  due  to  cooling  from 
the  fabricating  temperature.  He  considered  an 
idealized  composite  consisting  of  a  long 
sapphire  cylinder  surrounded  by  an  aluminum 
matrix.  This  is  equivalent  to  the  long -compo¬ 
site-cylinder  model  used  by  Dvorak  et  al.  This 
model  implies  that  the  thermal  strain  along 
the  cylinder  axis  is  a  constant,  i.e.  that 

dl 

e.  =  "  =  constant 
/ 

where  t  is  the  length  of  the  cylinder.  This  sim¬ 
plifies  the  procedure  of  obtaining  the  radial 
md  the  tangential  strains  which  was  done 
using  incompressibility  and  boundary  condi¬ 
tions.  The  effective  strain  e  at  the  Al-a-Al203 
interface  was  found  to  be  1.6%. 

However,  the  above-mentioned  composite 
models  do  not  give  an  accurate  description 
of  the  plastic  strain  state  in  the  short-compo¬ 
site-cylinder  model.  A  short-cylinder  model 
nearly  duplicates  the  situation  in  the  whisker 
and  platelet  Al-SiC  composites  at  present 
produced. 

The  purpose  of  this  work  was  to  determine 
experimentally  the  magnitude  of  the  local 
plastic  strain  produced  in  the  aluminum 
matrix  around  a  short  SiC  cylinder  during  a 
thermocycle  and  also  to  estimate  the  extent 
of  the  plastic  zone  around  the  cylinder.  Also 
an  effort  was  made  to  develop  a  theoretical 


model  of  the  plastic  zone  around  the  short 
cylinder. 


2.  EXPERIMENTAL  PROCEDURE 

The  small  particle  interspacing  (several 
micrometers)  in  a  commercially  available 
Al-SiC  composite  makes  it  impossible  to 
measure  the  local  plastic  strain  at  the  Al-SiC 
interface  directly .  Thus  a  composite  model 
consisting  of  an  SiC  cylinder  embedded  in  an 
aluminum  matrix  was  fabricated  so  that  direct 
strain  measurements  could  be  attempted. 

Aluminum  of  99.99%  purity  (to  eliminate 
the  influence  of  the  alloying  elements)  and 
commercial  carborundum  were  used  to  pro¬ 
duce  the  composite  model.  Selected  properties 
of  these  materials  are  given  in  Table  1 . 

Platelets  of  SiC  were  separated  from 
carborundum  conglomerates  that  are  used  in 
the  production  of  abrasives.  These  platelets 
were  spark  planed  on  an  electric  discharge 
machine  to  plates  approximately  1  mm  thick. 
These  flat  plates  were  cut  into  rectangular 
rods  approximately  1  mm  X 1  mm..  After  this, 
each  rod  was  spark  machined  to  cylinders  of 
about  1  mm  diameter. 

Pure  aluminum  rods,  12.5  mm  in  diameter 
in  the  as-received  condition,  wi  ■>  cut  into 
studs  37  mm  long.  Two  alumina.!  studs  and 
one  SiC  rod  were  assembled  together  and  put 
in  a  specially  built  compaction  die,  where 
they  were  hot  pressed  to  produce  one  com¬ 
pact.  Compaction  was  done  on  the  Instron 
testing  machine.  During  the  entire  compaction 
cycle  (Table  2),  a  vacuum  of  about  10~^  Torr 
was  maintained  using  a  mechanical  vacuum 
pump. 


TABLE  1 

Selected  properties  of  the  aluminum  and  SiC  used  in  this  paper 


.Material 

Yield 

strength 

(MPa) 

Elastic 

modulus 

(MPa) 

Hulk 

modulus 

(MPa) 

I’oisson 's  Thermal  c.xpansion 
ratio  coefficient 

(m  m“'  K”* ) 
in  the  temperature 
range  293-  773  K 

Melting 

point 

(K) 

Reference 

Al 

SiC 

ll,?,! 

.3  1..!  ■  (len.sion) 

62  X  10^ 
483  X  10^ 

57.5X103 
96.6  X  103 

0.31  28x10“*^ 

0.19  3xl0‘‘’ 

933 

3373 

110,  11] 
(12,  13) 

1.3^0  * 


(compression ) 


TABLE  2 


duced  where  Jcpss  is  equivalent  to  the  product 
Flow  chart  of  the  compaction  of  an  Al-SiC  composite  of  slip  band  density  N  and  the  amount  of  slip 

cylinder  S,i.e. 


Heating 

Heating  time 

=  1  h 

Temperature 

*  843  K 

Compression 

Cross-head  speed 

0.1  cm  min' 

Total  travel 

10  mm 

Holding  at  constant  load 

Time 

2  h 

Pressure 

2.28  MPa 

Incremental  cooling  under  pressure 

Time 

12  h 

Pressure 

2.28  MPa 

Fig.  1.  A  schematic  diagram  of  the  fabrication 
sequence  of  the  Al-SiC  composite  model. 


The  central  portion  of  the  compacted 
sample  was  sliced  in  the  transverse  direction 
into  disks  1  mm  thick  using  the  electric 
discharge  machine,  set  at  a  low  power,  and 
each  disk  contained  an  SiC  cylinder  very  close 
to  the  perfect  center  (Fig.  1).  All  Al-SiC 
disks  were  metallographically  and  then 
electrolytically  polished  to  remove  the  thin 
cold -worked  surface  layer  of  aluminum. 

The  method  adopted  for  the  evaluation  of 
plastic  deformation  was  based  on  the  direct 
observation  of  slip  bands  on  the  polished 
surface  of  the  sample  around  the  SiC  particles. 
The  amount  of  slip  is  a  characteristic  of  the 
amount  of  plastic  deformation  (when  defor¬ 
mation  occurs  by  slip)  in  a  crystalline  solid. 
The  plastic  strain  can  be  evaluated  if  the 
number  of  slip  bands  and  the  displacement  on 
each  band  are  known.  The  concept  of  com¬ 
bined  plastic  shear  strain  t^pss  has  been  intro- 


y,^,,  =  KNS  (1) 

where  K  is  a  coefficient  which  takes  into 
consideration  different  crystallographic  situa¬ 
tions.  A  detailed  treatment  of  is  given  in 
Appendix  A.  Thus,  the  method  reduces  the 
data  collection  to  the  measurements  of  slip 
band  densities  and  their  heights  in  the  area  of 
interest. 

The  electropolished  Al-SiC  disks  were 
separated  into  three  groups;  A,  B  and  C.  Each 
group  was  heated  to  about  823  K  and  then 
cooled  as  follows:  group  A,  furnace  cooled; 
group  B,  air  cooled;  group  C,  quenched  in 
alcohol. 

Since  the  surface  of  each  disk  had  a  high 
quality  polish,  slip  bands  could  be  observed 
around  the  SiC  in  an  optical  microscope. 

Slip  band  density  and  height  measurements 
were  obtained  using  a  Zeiss  interference 
microscope.  Areas  containing  slip  bands  were 
photographed  in  white  light  and  in  green 
monochromatic  light.  Pictures  taken  in  white 
light  gave  the  actual  image  of  the  slip  bands. 
Pictures  of  the  same  areas  taken  in  mono¬ 
chromatic  light  gave  interference  fringe 
patterns  (Figs.  2-4).  Thus,  the  correlation 
between  slip  bands  and  interference  fringes 
was  established. 

When  a  furrow  is  present  in  the  plane 
surface  of  the  object,  the  straight  interference 
bands  are  deflected  by  the  furrow  and  the 
furrow  depth  t  can  be  determined  from  the 
deflection: 


where  d  is  the  deflection  of  the  interference 
band,  b  is  the  band  interval,  taken  as  the 
distance  from  the  middle  of  one  band  to  the 
middle  of  the  next  band  and  X/2  =  0.27  pm 
for  the  thallium  light.  Band  deflection  can  be 
estimated  to  be  one-tenth  of  the  band  interval. 
Thus  the  height  measurements  can  be  as 
accurate  as  ±  0.1  X  0.27  =  ±  0  027  pm. 

The  slip  band  density  was  determined  by 
using  the  “mesh”  method  which  consisted  of 
the  following.  Pictures  of  slip  band  images 
Emd  interference  patterns  were  enlarged  to 
10  cmX  12.5  cm  on  high  contrast  photo¬ 
graphic  paper.  After  that,  a  transparent  plastic 
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TABLE  3 

Combined  plastic  shear  strain  7cpss  g''°up  A  samples  (furnace  cooled) 


_• _ _ 

<!LI'CT 

' 'A  ' 


Sample 

Path^  Combined  plastic  shear  strain^ 

7cpss  (^e)  /or  the  following  distances  from  the  Al-SiC  interface’^ 

6.25  X  10  ~ mm 

18.75  X 

'~mm  31.25X10  ~mm 

43. 75  X  10  ~mm 

56.25  X  10  ~  mm 

y,' 

’TT 

s'. 

1 

1 

1 

2 

0.86 

0.61 

0.35 

0.40 

0.18 

0.22 

0.09 

0.06 

0.026 

0.07 

I7 

4 

k  • 

1 

3 

0.41 

0.140 

0.24 

0.08 

1 

4 

0.26 

0.21 

0.19 

0.16 

adi 

1 

5 

0.42 

0.32 

0.32 

0.13 

0.03 

• 

1 

1 

6 

0.51 

0.41 

0.30 

0.08 

m' 

-1 

2 

1 

0.30 

0.21 

0.06 

2 

2 

0.45 

0.28 

3 

1 

0.43 

0.24 

0.22 

0.22 

0.22 

i'-'f-Z- 

■ 

3 

2 

0.26 

0.27 

0.22 

0.11 

A  a  JL  jj 

3 

3 

0.39 

0.22 

0.26 

0.16 

3 

4 

0.46 

0.18 

0.23 

X  ' 

3 

5 

0.45 

0.26 

1  . ' - ' 

V 

•a' 

4 

1 

0.43 

0.24 

0.04 

V  / 

‘«V-  !.•» 

CJh 

''iT) 

Average 

0.44 

0.24 

0.21 

0.16 

0.05 

■,!  - 

r-'.N 

4. 

“The  path  is  the  route  in  the  radial  direction  from  the  Al-SiC  interface  into  the  matrix,  which  was  selected  for 
slip  band  density  and  height  measurements.  Selection  was  based  on  the  amount  of  slip  that  occurred,  and  those 
routes  were  selected  in  which  the  amount  of  slip  bands  appeared  to  be  the  greatest. 

^The  given  strain  is  actually  the  result  of  the  product  NS.  It  does  not  include  the  coefficient  K  which  was  found 
to  be  equal  to  3  (see  Appendix  A).  Values  of  7cpss  that  incorporate  K  are  given  in  Table  6. 

'The  distance  given  here  is  the  distance  from  the  Al-SiC  interface  to  the  centers  of  the  first,  second  etc.  squares. 
Therefore,  values  of  Jcpss  represent  the  strain  over  the  entire  square. 

TABLE  4 

Combined  plastic  shear  strain  7cpss  '*1  group  B  samples  (air  cooled) 


Sample 

Path^ 

Combined  plastic  shear  strain^ 

Tepss  (^)  for  the  following  distances  from  the  Al-SiC  interface’^ 

6.25  X 

18.75  X 

31.25  X 

43.  75  X 

56.25  X 

68.  75  X 

:  „  ,1 

• 

10  -’mm 

10  ^  mm 

lO'^mm 

10  ^  mm 

10  ~  mm 

10  -  mm 

1 

1 

0.34 

0.32 

0.35 

0,13 

1 

2 

0.25 

0.25 

0.17 

0.19 

0.13 

0.11 

1 

3 

0.56 

0.31 

0.19 

0.11 

1 

4 

0.47 

0.55 

i.s'V/'Zi 

1 

5 

0.25 

0.07 

0.014 

0.08 

2 

1 

0.23 

0.17 

5.' 

■1 

1 

0.65 

0.41 

1 

2 

0.47 

0.22 

•1 

3 

0.54 

0.29 

0.1 

4 

1 

1 

r> 

0.31 

0.48 

0.21 

0.25 

0.25 

9. 

Avora^o 

0.41 

0.27 

0.18 

0,12 

0.13 

0.11 

■*  The  path  is  the  route  in  the  radial  direction  from  the  Al-SiC  interface  into  the  matrix,  which  was  selected  for 
slip  band  density  and  height  measurements.  Selection  was  based  on  the  amount  of  slip  that  occurred,  and  those 
routes  were  selected  in  which  the  amount  of  slip  bands  appeared  to  be  the  greatest. 

The  given  strain  is  actually  the  result  of  the  product  NS.  It  does  not  include  the  coefficient  K  which  was  found 
to  be  equal  to  3  (see  Appendix  A).  Values  of  7ci>ss  ^-^at  incorporate  K  are  given  in  Table  6. 

'The  distance  given  here  is  the  distance  from  the  Al-SiC  interface  to  the  centers  of  the  first,  second  etc.  squares. 
Therefore,  values  of  represent  the  strain  over  the  entire  square. 
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TABLE  6 

Values  of  7epss  incorporating  K  for  group  A,  B  and  C  samples 

Group  Combined  plastic  shear  strain  7cpss  (”  )  following  distances  from  the  Al-SiC  interface 


6.25  X 

18.75  X 

31.25  X 

43. 75  X 

56.25  X 

68.  75  X 

1 0~-  mm 

1 0~-  mm 

10~- mm 

10~  -  mm 

1  ()~-  mm 

10~~  mm 

A 

1.32 

0.72 

0.63 

0.48 

0.15 

B 

1.23 

0.81 

0.54 

0.36 

0.39 

0.33 

C 

0.99 

0.60 

0.3 

determination  can  be  evaluated  as 
±  (0.027/0.135)  X  100  ==  20%  where  0.135  pm 
is  taken  as  the  average  slip  band  height.  This  is 
obviously  a  large  error  range. 

3.2.  Hot-stage  observations 

3.2.1.  Heating 

Slip  bands  began  to  appear  even  before  the 
temperature  reached  373  K  and  were  forming 
as  widely  spaced  deep  lines.  As  the  tempera¬ 
ture  increased,  the  density  of  the  slip  bands 
increased  also,  reaching  an  apparent  maxi¬ 
mum  at  around  573  K.  At  around  533  K  a 
new  group  of  broken  lines  began  to  form  with 
no  relation  to  the  previously  formed  slip 
bands.  The  development  of  these  lines  became 
more  intensive  at  higher  temperatures.  When 
the  temperature  approached  about  the  653  K 
mark,  the  separation  between  the  matrix  and 
the  SiC  became  fairly  visible.  At  approxi¬ 
mately  693  K,  slip  bands  began  to  fade  and 
disappeared  almost  totally  by  the  time  the 
temperature  reached  853-873  K.  At  the 
maximum  temperatures,  the  surface  did  not 
appear  as  flat  as  it  did  at  the  beginning. 

3.2.2.  Cooling 

On  cooling,  previously  formed  slip  band 
patterns  began  to  show  up  again  but  they  did 
not  reat:h  the  size  and  extent  of  the  former 
slip  band  patterns.  I'he  development  of  the 
slip  bands  in  new  areas  was  not  observed. 

The  broken  line  which  encircled  the  entire 
SiC  particle  began  to  develop  radial  spokes 
that  connected  this  line  to  the  interface. 

When  the  sample  cooled  to  room  temperature, 
the  surface  of  the  sample  remained  rippled. 

Subsequent  thermocycles  did  not  show  any 
changes  in  slip  band  morphology.  A  broken 
line  around  the  SiC  cylinder  became  much 
more  clearly  defined  with  repeated  thermo¬ 
cycles.  The  matrix  appeared  to  be  separated 
from  the  particle  all  the  time. 


4.  DISCUSSION 

Composite  models  discussed  in  the  litera¬ 
ture  for  the  theoretical  determination  of 
stresses  and  strains  in  the  matrix  (spherical 
particle  in  a  matrix,  long  composite  cylinder 
etc.)  cannot  be  applied  to  a  short  composite 
disk.  Therefore,  an  attempt  was  made  to 
evaluate  the  radius  of  the  plastic  zone  and  the 
plastic  strains  around  the  SiC  particle  using 
the  approach  described  in  Appendix  B. 

The  resultant  profile  of  the  effective  plastic 
strain  e  around  the  SiC  particle  and  the  extent 
of  the  plastic  zone  are  shown  in  Fig.  5.  As 
can  be  seen,  p  ^  1.15  mm,  and  e  =  2,5%  at 
the  interface.  The  theoretical  plastic  strain  is 
higher  than  the  observed  plastic  strain,  and 
the  theoretical  plastic  zone  radius  is  smaller 
than  the  observed  plastic  zone  radius.  The 
theoretical  plastic  zone  size  is  estimated  on 
the  assumption  that  deformation  is  homogen¬ 
eous,  i.e.  uniform  around  the  SiC  particle.  In 
the  real  case  we  have  “bursts”  of  plastic  flow 
in  accordance  with  a  particular  crystallographic 
situation.  Thus,  the  extent  of  the  burst  of 
plastic  deformation  can  be  larger  than  that  of 
homogeneous  deformation  since  the  same 
amount  of  metal  flow  must  be  accommoda¬ 
ted.  The  explanation  of  the  difference  between 
observed  and  calculated  strains  can  be  given  as 
follows.  On  heating,  because  of  the  differen¬ 
tial  expansion  between  the  aluminum  and  the 
SiC,  the  aluminum  matrix  tends  to  pull  away 
from  the  SiC.  If  the  bonding  between  the 
aluminum  and  the  SiC  is  sufficient  to  resist 
the  pulling  action,  the  matrix  undergoes  de¬ 
formation  proportional  to  AT  of  the  cycle.  If, 
however,  the  bonding  is  weak,  the  matrix 
breaks  away  from  the  SiC  at  a  certain  temper¬ 
ature  T/, ,  and  subsequent  temperature  in¬ 
crease  has  no  effect  on  matrix  deformation 
since  it  is  now  free  to  expand. 
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Hot-Stage  observations  showed  that  slip 
band  formation  started  at  a  temperature  of 
less  than  373  K.  Accurate  evaluation  of  the 
temperature  at  which  slip  became  visible  was 
not  possible  because  the  thermocouple  in  the 
hot  stage  was  not  sensitive  enough  in  the 
temperature  range  between  298  and  373  K. 

It  was  notx  ed  that  more  slip  bands  were 
formed  on  the  heating  half  of  the  thermocycle 
than  on  the  cooling  half.  The  slip  band 
arrangement  on  the  cooling  half  ’■epeated  the 
arrangement  developed  during  the  heating 
half.  When  the  temperature  approached  about 
773  K.  slip  bands  began  to  disappear.  This  can 
be  e.xplained  from  the  surface  tension  point 
of  view. 

■Vnother  possible  explanation  of  the 
stoppage  of  slip  band  development  is  that  the 
aluminum  matrix  breaks  away  from  the  SiC 
and  continues  to  expand  freely  without  any 
restraint  from  the  SiC.  It  is  also  possible  that 
very  fine  slip  bands  still  continued  to  form 
(high  temperature  creep,  for  example),  but 
we  did  not  see  this  because  of  the  limitation 
in  resolution  of  our  optical  system.  It  should 
be  mentioned,  however,  that  the  disappearance 
of  the  slip  bands  corresponds  to  a  similar 
observation  made  by  Lammers  et  al.  [14]  in 
their  m  nitu  transmission  electron  microscopy 
investigation  of  Al-SiC  composites.  They 
observed  “slip  lines”  in  thermally  cycled 
transmission  electron  microscopy  foils,  and 
these  slip  lines  disappetired  at  high  tempera¬ 
tures. 

The  fact  that  slip  band  patterns  formed  on 
cooling  repeat  themselves  shows  that  the  same 
slip  systems  are  engaged  in  the  “reversed” 
deformation,  proving  at  least  partially  that 
there  is  a  certain  reversibility  of  the  plastic 
deformation  that  occurs  when  the  load  is 
reversed. 

The  broken  lines  that  form  around  SiC 
appeared  to  be  the  boundaries  of  recrystal¬ 
lized  grains.  If  this  is  indeed  the  case,  re¬ 
crystallization  took  place  at  the  interface  and 
occurred  very  rapidly  (2  3  min). 


.■>  CONCLUSION 

From  the  experimental  data  and  the  theo¬ 
retical  model  the  following  conclusions  were 
obtained. 


(1)  A  new  experimental  technique  for 
measuring  local  plastic  strains  was  developed 
and  utilized  in  the  /\l-^SiC  system. 

(2)  The  theoretical  treatment  of  the  short 
composite  cylinder  gave  the  distribution  and 
extent  of  the  plastic  strains,  which  are  in  fair 
agreement  with  the  experimental  results. 
There  is  a  disagreement,  however,  between 
the  experimental  and  theoretical  values  of 
plastic  strains  immediately  at  the  Al-SiC 
interface . 

(3)  An  Al-SiC  bond  is  a  very  important 
factor  influencing  the  plastic  deformation 
around  SiC  particles. 
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APPENDIX  A 

A.l.  Combined  plastic  shear  strain 

The  relation  between  the  slip  band  density 
N,  the  amounts  of  slip  or  displacement  on 
the  band  and  the  plastic  strain  y  has  been 
discussed  in  a  number  of  publications  [  Al- 
A3|.  In  their  classical  works,  Yamaguchi  et 
al.  I  Al  I  and  Brown  [A2]  found  direct  pro¬ 
portionality  between  the  strain  and  amount 
of  sli{>  during  plastic  deformation.  Cottrell 


[A3]  summarized  their  results  and  arrived  at 
the  expression  7  =  NS  for  the  case  when  there 
is  only  one  slip  system  in  operation. 

For  a  multiple-slip-system  case,  however, 
this  expression  must  be  modified  to  take  into 
account  the  different  crystallographic  orienta¬ 
tions  of  the  slip  systems. 

The  attempt  to  carry  out  this  modification 
is  offered  in  the  present  work  and  consists  of 
the  followdng:  the  highlights  of  the  one-slip- 
system  case;  a  description  of  the  rigorous 
approach  for  multislip  cases;  a  simplified 
approximation  with  the  introduction  of  the 
combined  plastic  shear  strain  7cpss  concept. 

Let  us  assume  first  that  slip  occurs  in  one 
slip  system  only.  Let  us  consider  an  imaginary 
block  (Fig.  Al)  that  is  cut  out  of  material  in 
such  a  way  that  slip  planes  are  perpendicular 
to  the  facets  of  the  block.  This  block  is 
oriented  in  such  a  way  that  its  height  is 
measured  along  the  z  axis. 

Let  us  suppose  that  a  shear  stress  is 
applied  to  the  block  which  causes  the  block 
to  slip  as  shown  in  Fig.  Al . 

By  definition,  the  shear  strain  in  this  case 
is  expressed  as 


where  S,  is  the  slip  or  the  displacement  in  the 
ith  band*.  (Here  S,  is  analogous  to  the  furrow 
depth  t  in  eqn.  (2).)  If  n  is  the  total  number 
of  slip  bands  involved  in  the  slip  and  Sj,  is 
the  average  value  of  slip  per  one  band,  then 

AS,  =  X  S- 

1  =  1 


As  can  be  seen  from  Fig.  Al,  n  is  the  number 
of  intersects  that  the  2  axis  makes  with  the 
slip  bands  and  is  equal  to  the  total  number  of 
slip  bands  along  the  2  direction.  If  the  number 
N  of  slip  bands  per  unit  length  along  z  is 
known  (i.e.  the  slip  band  density),  then 


where  N  (cm'^)  is  the  slip  band  density  along 
the  2  direction.  Now  the  total  strain  7,^  can 
be  expressed  [A3]  as 

2 


^ 


»  i 


where  AS,  is  the  total  displacement  in  the 
y  direction  and  Z  is  the  height  of  the  block. 

The  total  displacement  AS,  is  equEil  to  the 
sum  of  the  displacements  on  each  slip  pleme, 
i.e. 

AS,  =  i  S, 

i=l 


=  NSy 

N  can  be  evaluated  as  the  total  slip  band 
length  divided  by  the  total  area.  To  do  this 
we  can  take  eqn.  (Al)  and  multiply  the 
numerator  and  the  denominator  by  the  width 
L  of  the  block: 


^  aw 

Of 

3 


/ 

7' 

"  _  i 

L'Si  7 

i 

/ 

T- 


Fig.  Al,  Hypothetical  block  with  one  slip  system. 
Sj  is  the  displacement  on  the  slip  plane;  AS,  is  the 
total  displacement  along  the  y  axis. 


_  nSy  L 
z  L 

^  {nL)Sy 
zL 

‘Generally,  the  slip  direction  does  not  have  to  be 
perpendicular  to  the  .v  plane  as  shown  in  Fig.  Al  . 
However,  despite  this  and  other  rather  crude  assump¬ 
tions  that  are  made  later,  the  purpose  is  to  show  the 
complexity  of  the  rigorous  treatment. 


—I 


‘  O  •L  ■  •». 
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Fig.  A2.  Hypothetical  block  with  two  operational  slip 
systems  1  and  2. 


Here,  nL  is  the  total  length  of  the  slip  bands 
on  the  side  of  the  block,  zL  is  the  total  area 
of  this  side  and  Np^  (cm"^)  =  nLjzL  is  the  slip 
band  density.  Thus, 

7>4  =  NSy 

In  the  case  when  more  than  one  slip  system  is 
involved  in  plastic  deformation,  slip  bands 
are  not  parallel  but  instead  each  slip  system 
has  its  own  orientation*  (Fi^;.  A2).  In  this 
situation  the  equation  for  shear  streiin  has  to 
be  modified  to  take  into  account  the  different 
orientations  of  the  slip  bands.  The  rigorous 
way  to  do  this  would  be  as  follows. 

(1)  Determine  the  plastic  strain  7^  ^  using 
the  same  approach  described  previously  for 
slip  system  1. 

(2)  Rotate  the  coordinate  system  ,  yi , 
to  align  it  with  system  ;c2,  y2>  ■^2  where  22  is 
the  direction  perpendicular  to  slip  bands  in 
slip  system  2.  In  the  coordinate  system  X2,y2, 
22,  determine  7^,  ,,  using  the  relation 

yy:^:  =  (A2) 

where  is  the  rotation  matrix  and 

jjj  is  the  strain  tensor  in  the  aci  ,  yi ,  2i  coordi¬ 
nate  system  and  is  given  by 


0  0 


*Here,  Ihe  hypothetical  assumption  is  made  that 
slip  planes  from  the  second  slip  system  are  also 
perpendicular  to  the  facets  of  the  block. 


(3)  Combine  the  terms  determined  in  step 
(3)  with  the  terms  from  step  (1)  having  the 
same  indices.  Now  we  have  obtained  the  total 
strain  tensor  e,y  in  the  Xi,yi,Zi  coordinate 
system  due  to  all  slip  systems: 


(4)  Determine  the  principal  components  of 
the  strain  tensor  e,y  solving  the  following 
equation; 

r-X  0  0  1 


for  X.  The  three  roots  Xj,  X2  and  X3  are  the 
three  principal  strains  ,  €2  and  C3  respec¬ 
tively. 

(5)  Determine  the  effective  strain 
21/2 

e  =  — {(€1-62)^  +  (ea  -  63 +  (63-61)=*}^=^ 

As  a  result  of  all  these  steps  we  finally  get 
one  data  point  on  the  plot  e  =  f(r).  This 
procedure  is  simple  in  principle  but  compli¬ 
cated  to  implement,  since  the  crystallographic 
orientation  has  to  be  determined  for  each 
grain.  As  an  alternative  way  of  estimating  the 
amount  of  plastic  deformation  that  occurred 
around  the  SiC  particle,  the  following  approxi¬ 
mation  is  offered. 

First,  the  concept  of  combined  plastic 
shear  strain  is  introduced.  We  shall  define 
Tcpss  as 

Tcpss  =  KNS 

where  S  is  the  average  step  height.  In  general, 
the  slip  plane  is  not  perpendicular  to  the 
surface  and  therefore  the  coefficient  K  is 
incorporated.  K  is  a  coefficient  which  will  be 
evaluated  later.  N  is  the  slip  band  density, 
which  is  obtained  from  the  total  slip  band 
length  of  all  slip  systems  divided  by  the  area 
where  the  total  slip  band  length  is  the  length 
of  all  the  slip  bands  enclosed  in  a  selected 
area  {e.g.  a  square  mesh  of  some  net).  The 
selection  of  the  size  of  the  area  is  based  on 
considerations  of  the  scale.  It  is  realized  that 
the  combined  plastic  shear  strain  7cpss  has  no 
direct  correlation  with  the  effective  strain  e, 
the  octahedral  shear  strain  7„ci  or  any  other 
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specific  strain.  However,  the  relative  simplicity 
in  the  determination  of  y^pss  and  the  fact  that 
Tcpss  takes  into  account  direct  proportionality 
between  the  amount  of  plastic  deformation 
and  the  slip  band  density  makes  it  very  attrac¬ 
tive  for  use  as  a  criterion  for  plastic  deforma¬ 
tion  in  our  case  and  in  other  similar  cases  in 
general. 

A. 2.  Determination  of  the  coefficient  K 
The  coefficient  K  takes  into  account  the 
different  crystallographic  orientations  of  slip 
systems.  Let  us  consider  the  case  when  we 
have  three  slip  systems  operating,  i.e.  we  see 
three  slip  band  groups.  The  individual  slip 
systems  contribute  the  strain  ATjS, ,  N2S2  and 
N3S3  to  the  total  i.e. 

Tcpss  ~  ^3^3 

where  Alj,  Nj  andN3  are  the  slip  band  densities 
and  Si,  S2  and  S3  are  the  average  displace¬ 
ments  on  the  band  within  each  slip  system, 
given  by 


■Smi 
sin  0i 


and 


sin  O3 

where  ,  iSm2  .§^3  are  the  average 
heights  of  the  slip  bands  in  the  three  slip 
systems,  measured  perpendiculeirly  to  the 
surface,  and  0i ,  O2  and  O3  are  the  angles 
between  the  slip  planes  and  the  surface 
plane.  Thus, 


X 

Tcpss  ~  ^iSmi  .  „ 

sm  0i 


+  A^3Sm3 


1 

sin  O3 


Now  we  shall  assume  that  the  amounts  of 
strain  NiS^i ,  N2SM2  approxi¬ 

mately  equal  because  of  the  symmetry  of  the 
specimen  and  provided  that  no  voids  or  cracks 
are  created.  Furthermore, 


1  1 

- - .  - - 

sin  di  sin  O2 


or 

Tcpss 

where 


111 

—  —  -f.  -  —  4- - 

sin  dj  sin  O2  sin  O3 


Let  us  evaluate  the  range  within  which  K  can 
change  by  considering  a  unit  triangle  (Fig. 
A3).  Let  the  directions  shown  on  the  triangle 
correspond  to  the  poles  of  the  surface  of  the 
sample.  Then  the  values  of  K  would  be  as 
shown  in  Table  A1 .  The  average  value  of  K 
is  about  4.5  when  three  slip  systems  are 
operating.  Thus,  the  average  correction  factor 
per  one  slip  system  is  4.5/3  =  1.5.  In  general, 
slip  band  densities  were  measured  along 
directions  going  through  areas  where  two  slip 
systems,  j.e.  two  sets  of  slip  bands,  were 
present.  Thus  for,  two  slip  systems, K  =  2X 
1.5  =  3. 


Fig.  A3.  Some  orientations  of  the  sample  for  which 
values  of  K  were  determined. 


TABLE  A1 

Values  of  K  corresponding  to  the  orientations  shown  in  F.,;.  A3 


Surface  plane  pole  (111)  (110)  (100)  (210)  (211)  (221)  (310)  (311)  (320)  (321) 

K  for  ( 11 1 }  slip  plane  3.2  3.7  3.7  3.66  .5.1  5.9  3.6  4.2  3.1  4.9 
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APPENDIX  B 

Among  the  several  composite  models  des¬ 
cribed  in  the  literature,  there  is  one  that 
appears  to  be  close  to  our  case;  this  is  known 
as  the  long-composite-cylinder  model.  It 
consists  of  long  continuous  fiber  surrounded 
by  the  metal  matrix.  This  model  implies  that 
plastic  strain  in  the  direction  of  the  fiber  (the 
z  direction)  resulted  from  unequal  expansion 
or  contraction  of  fiber  and  that  the  matrix 
is  a  constant,  i.e.  e,  =  dl/l  =  constant  where  / 
is  the  length  of  the  fiber. 

Therefore,  the  total  plastic  strain  state 
becomes  a  plane  strain  case  which  enables  two 
other  plastic  strains  e,.  and  eg  to  be  found 
without  difficulty.  Although  short  composite 
cylinders  have  the  same  symmetry  as  long 
composite  cylinders,  the  strain  in  the  z 
direction  is  not  a  constant  but  is  in  general  a 
function  of  two  variables:  (1)  the  distance 
from  the  interface  and  (2)  the  position  along 
the  z  axis.  Since  +  constant,  the  long- 
cylinder  approach  cannot  be  used  in  our  case. 

VVe  shall  discuss  the  case  shown  in  Fig.  Bl. 
Let  us  consider  an  Al-SiC  disk  that  cools 
from  a  temperature  T.  As  a  result  of  cooling, 
a  certain  amount  of  aluminum  will  be  pushed 
back  because  of  the  differential  shrinkage  of 
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SiC  and  aluminum.  If  we  did  not  have  SiC  in 
the  center  of  the  aluminum  disk,  the  alumi¬ 
num  ring  would  shrink  without  any  restraint 
and  the  size  of  the  bore  would  be  reduced  by 

U„  =  Aa  -  a  Aa  AT 

where  a  is  the  radius  of  the  SiC,  Aa  (=5=  25  X 
10'®  K“^)  is  the  difference  between  the  co¬ 
efficient  of  thermal  expansion  of  aluminum 
and  that  of  SiC,  AT  500  K)  is  the  tempera¬ 
ture  interval  and  [4  is  the  displacement  at  the 
interface.  The  reduction  in  the  radius  of  the 
bore  does  not  occur  when  SiC  is  present  in 
the  center  of  the  aluminum  disk.  The  amount 
of  aluminum  that  is  not  allowed  to  go  towards 
the  center  of  the  bore  would  cause  plastic 
flow  of  the  adjacent  matrix  in  all  directions 
away  from  the  SiC.  Thus,  some  of  the  matrix 
will  spill  out  (as  shown  in  Fig.  Bl).  As  we  go 
further  away  from  the  interface,  the  plastic 
flow  is  restricted  by  the  matrix  that  surrounds 
the  central  portion .  Two  regions  can  be 
considered:  the  plastic  region  that  we  assume 
is  adjacent  to  the  SiC  and  the  elastic  region 
that  encloses  the  plastic  region.  Let  us  make 
the  following  assumptions. 

(1)  In  the  plastic  zone  the  matrix  is  a  per¬ 
fectly  plastic  material,  i.e.  no  work  hardening 
occurs. 

(2)  At  the  starting  temperature  T  of  773  K, 
SiC  and  aluminum  are  just  in  contact  with 
one  another. 

(3)  The  profile  of  the  aluminum  which  has 
spilt  out  is  a  straight  line. 


Fig.  Bl .  Schematic  representation  of  the  plastic-elastic  shell  around  the  SiC. 


(4)  The  plastic -elastic  front  is  also  a 
straight  line. 

Now  let  us  consider  the  radius  p  of  the 
plastic  zone.  It  is  a  function  of  the  vertical 
position  of  the  transverse  plane  in  which  p 
is  considered.  According  to  our  assumption, 

Pq-Qz 

(the  equation  of  the  straight  line)  where  p^  is 
the  plastic  zone  radius  in  the  “0”  plane, 
which  is  the  plane  of  symmetry  between  the 
upper  and  lower  halves  of  the  composite 
disk.  To  determine  q  we  let  z  =  \  H,  i.e.  the 
top  of  the  “H”  plane  is  considered.  When 
z  =  \H,  p  =  p„.  Therefore 


From  this 


q  =  2(Po  -  Ph  ] 


and  now 


Po  2(Po 


Let  us  consider  now  the  portion  of  the  matrix 
that  spills  out.  From  the  triangle  shown  in 
Fig.  B2  we  can  write 


U,  =  Ah 


Ph  -r 
Pi  I  -a 


where  f4  is  the  vertical  displacement.  Gener¬ 
ally  speaking,  the  vertical  displacement  is 
a  function  of  z  and  r,  i.e. 

Uz  =  <t>{f{z),g{r)] 

Equation  (B2)  gives  only  an  expression  for 
when  the  “H”  plane  is  considered.  The 
problem  thus  is  to  find  a  general  expression 
for  U^.  To  do  this,  let  us  replace  p„  in  eqn. 


H"  PLANE 


z  =  1/2  H 


Fig.  B2.  “Spill-out”  portion  of  the  aluminum  matrix. 
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(B2)  with  a  general  expression  for  p  (eqn. 
(Bl))  when  z  =  \H 

-2{po— Pii)\  —  r 
U,^  Ah  -  f 

p  2(pQ-p„)i-a 

Ah  =  \H  AaAT 


U,  AaAT 


Po  2(Po  P//I2  ^ 

Po  2(Po  “  P// )  2  ~ 


If  we  simplify  eqn.  (B3)  it  reduces  to  eqn. 
(B2),  but  let  us  leave  eqn.  (B3)  in  its  form  and 
examine  it  in  more  detail.  First  let  us  substi- 
tude  H  ~  1  mm  in  eqn.  (B3),  since  this  is  a 
convenient  way  of  simplifying  this  expression 
and  H  -  1  mm  happens  to  be  the  actual 
height  of  our  composite  disks.  Now  we  can 
rewrite  eqn.  (B3)  as 

1  Po  ^ 2(Po  P/zlo  ~r 

U,-=lAa  AT  --  -  --  (B4) 

Po  2(Po  P//)2 

for  z  =  |//  or  if  we  normalize  z'  -  zjH  =  \ . 
Let  us  check  the  values  of  [4  at  different  z 
values.  First,  let  2  =  |  (“H”  plane),  and  so 

Po~2(Po-  Ph)\-  a 
Then,  when  r  =  a, 
f4  =  I  Act  AT  =  Ah 
and,  when  r  =  pn, 
t4  =  0 

Let  2  =  0  (“0”  plane),  and  so  =  0,  since 
the  “0”  plane  is  the  plane  of  vertical  sym¬ 
metry  .  At  this  point  we  can  see  a  certain 
logic  in  eqn.  (B4),  which  enables  us  to  suggest 
the  general  expression  for  in  the  following 
form: 


U,=z'  AaAT 


Po-2(Po-p„)2'  -r 
Po-2(Po-P//)z'  -  a 


(since  /f  =  1).  It  is  realized  that  eqn.  (B5) 
lacks  rigorous  proof.  However,  since  it  in¬ 
corporates  our  boundary  conditions  and 
since  we  do  not  have  any  additional  informa¬ 
tion,  it  is  reasonable  to  adopt  eqn.  (B5)  and 
to  see  what  results  we  can  obtain  and  then  to 
judge  the  validity  of  this  expression.  Let  us 
consider  the  cylindrical  coordinate  system  in 
which  the  2  axis  coincides  with  the  2  axis  of 
our  sample.  By  definition 


-  .'-V-  A .  - . 


e 


To  find  the  constant  C  we  use  the  boundary 
condition  Ur=a  =  U^=  a  Aa  AT: 


M  C  N  „ 

U„=  a  Aa.  AT  = - oH - 1 — 

2  a  3 

M  ^  N  , 

C^aU„+~a^ - 

2  3 


=  Ac  + 

Po  2(po  Ph)^  0 


+  z Aa AT 


2{Po-PH)(a  -r) 


[po-2(po  -Pn)z-a] 
Incompressibility  requires  that 

Cr  +  e„  +  e^,  -  0 


dU,  U, 

+  _1  +  e,  -  0 

dr  r 

Let  us  consider  the  “0”  plane,  where  2  =  0; 
then 


e^.  =  Aa  AT 


Po-r 

Po-a 


Equation  (B6)  then  becomes 

dt/r  ^4  Po~f 

L  +  Aa  AT  =0 

dr  r  Po~ci 

dU,  U, 

— -  +  M  -  Nr  =  0 

dr  r 

where 


Now  eqn.  (BIO)  becomes 

M  U^a  Ma^  Na^  N  , 

17,  =  --r+— + - +-r2  (Bll) 

2  r  2r  3r  3 

Now  we  can  find  the  strains  e„  €y  and  for 
the  “0”  plane: 

dl7,  M  U^a  Ma^  Na^  2 

e  =  — ^  = - - - + - +  -  ATr 

^  dr  2  r2  2r2  3r2  3 

M  U^a  Ma^  Na^  N 

r  “  2  ^  r^  2r2  3r2  "s 


e,  =  Af -Nr 


(B12) 


Let  us  consider  the  tip  of  the  plastic -elastic 
interface,  i.e.  r  =  p^.  Then  (e^  =  0  from  eqn. 
(B7)) 


fr  =  -  ey 


(B13) 


where  Up^  is  the  displacement  of  the  tip  of  the 
plastic -elastic  front: 

=  t/p„ 

Each  point  of  the  plastic-elastic  front  is  in 
the  state  of  incipient  yielding,  i.e. 


M  ^  Aa  AT 


Po -a 


N  =  Aa  AT - 

Po^a 

The  solution  of  eqn.  (B9)  is 

M  C  N 

Ur  =  -  -  r  +  -  +  -r^ 

2  r  3 


(BIO) 


where  e  is  the  effective  strain  at  the  interface 
and  Cp  is  the  yielding  strain  taken  from  the 
yielding  condition  of  the  tensile  sample  in 
uniaxial  loading: 

Oy  -  Eeo 

Therefore 


*From  here  on  we  omit  the  superscript  prime  in  z' 
and  use  z  instead  to  prevent  confusion,  but  it  should 
be  borne  in  mind  that  the  value  z  is  normalized  with 
respect  to  //. 


11.73 
62X10^ 
^  2X10' 


where  Oy  =  11.73  MPa  and  E  -  62  X  10^  MPa 
for  pure  aluminum.  Thus  we  have  a  system 
of  equations; 


U. 


^r(r  =  p„'> 


^0(r  =  P„) 


P,i 

Po 


U, 


_ 

Po 
=  0 


(B14) 


Applying  the  von  Mises  yielding  condition, 
21/2 


{(Cl  +  (61-63)"  +  (6.3  -63)"}^" 


where  63, 62  and  63  are  the  principal  strains. 
Generally  speaking,  the  strains  e^,  6y  and  6,  do 
not  coincide  with  the  principal  strains  63 , 62 
and  63,  because  of  the  existence  of  the  shear 
component  in  the  total  strain  tensor,  which  is 
given  by 

dG.  dU, 

The  components  yor  and  y^Q  are  equal  to  zero 
because  of  the  symmetry  around  SiC. 

df4  _  _  AaATz  _ 

dr  Po -2(Po -PH)2-a 

for  the  “0”  plane;  when  2  =  0,  dU^/dr  —  0.  To 
ev2iluate  dU^ldz,  let  us  consider  the  region  of 
the  matrix  adjacent  to  the  “0”  plane  as  shown 
in  Fig.  B3.  Point  A  is  displaced  to  point  A’, 
and  point  B  to  point  B  .  (B  is  a  mirror  image 
of  A.)  The  displacements  C/^a  and  C4b 
equal  because  of  the  symmetry  and  from  the 
assumption  that  no  voids  are  created.  The 
displacement  varies  as  we  go  along  the  2 
axis  and  presents  a  continuous  function  of  2. 
Since  the  function  is  continuous  2md  sym- 


- "0”  PLANE 


Kiy  B.3.  Uisplacomonts  AA'  and  BB'  of  the  matrix 
.symmetrical  about  the  “0”  plane. 


metrical  about  the  “0”  plane,  it  goes  through 
its  local  maximum  or  minimum  point  when  it 
crosses  the  “0”  plane  and  therefore  dUjdz  =  0 
when  2  =  0.  Thus  we  found  that,  within  the 
“0”  plane,  all  the  shear  components  of  the 
strain  tensor  are  equal  to  zero,  which  means 
that  6r  =  63 , 6o  =  62  and  =  63  for  the 
“0”  plane.  Thus,  if  we  go  back  to  the  von 
Mises  condition  of  yield,  the  effective  strain 
6  can  now  be  expressed  as 
21/2 

e=  — {(6,-60)"  +  (6„- 6,  )"  +  (6,-6,)" 


and,  from  eqn.  (B14), 


3  IV 
=  1.155 
However, 


Po 

Po 


■J  KITF" 


e  =  60 
=  2X10''* 
U. 


=  1.155  ■ 
Po 


Therefore 


— ’  =  1.7  XIO"'* 
Po 


(B15) 


From  eqn.  (Bll)  we  get 


M  U„a  Ma^  Na^  N  „ 

-  Po  +  -+ - +-Po^ 

2  Po  2po  3po  3^° 


and,  substituting  for  M  and  N,  we  obtain 
AaArjPo"+apo+a"  (a+Po)Po 


Up.- 


Po 


+  a" 
(B16) 


Equations  (B15)  and  (B16)  can  be  solved  as  a 
system  of  two  equations  with  two  unknowns, 
(4,  and  Po.  Solving  them,  we  obtain  Po  =1.15 
mm.  In  order  to  find  P/,  we  can  follow  the 
same  procedure.  However,  for  the  “//”  pleme 
case  the  strains  6,,  6,,  and  6,  are  not  equal  to 
63, 62  and  63.  Therefore  the  procedure  has  to 
be  modified.  The  new  procedure  would  be  to 
find  an  expression  for 


dt/.  dt/, 

r  =  *  - 

dr  d2 


7.- 


f 
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for  2  I  and  then  to  solve  the  equation 


e,  -  X 

0 

7,-. 

0 

eg  -  A 

0 

7--r 

0 

e,  -  A 

This  is  a  cubic  equation.  Three  roots  give  the 
values  for  the  principal  strains; 

Xi  =  Cl 
^2  ~  ^2 
and 

^3  =  ^3 

This  is  simple  in  principle  but  very  complex 
to  implement  for  the  follovving  reason.  The 
individual  components  of  this  tensor  are  of 
the  type  expressed  by  eqns.  (Bll)  and  (B15). 
If  the  coefficients  of  the  cubic  equation  were 
numerical,  then  we  could  solve  it  using  the 
trial-£md-error  method.  In  our  case,  all  the 
coefficients  have  a  general  expression .  The 
general  solution  of  the  cubic  equation  has  the 
form 

r  n  I  /n\3ll/2ll/3 


3il/2-il/3 


Each  q  and  p  would  be  a  combination  of  the 
expressions  for  e„  eg,  and  Because  of 
the  difficulty  of  determining  in  a  rigorous 
way  we  can  use  an  approximate  solution  by 
treating  strains  £„  £g  and  e.  as  the  principal 
strains  e'l ,  62  and  63.  This  implies  that  we 
ignore  the  influence  of  the  shear  strain  7-,.. 

It  is  realized  that  the  resulting  value  for  p„ 
will  not  be  a  true  value  but  it  will  at  least  give 


an  idea  of  the  extent  of  the  plastic  zone  on 
the  surface  of  the  sample.  Additionally,  if 
the  magnitude  of  7^^  is  small,  the  deviation  of 
e„  £()  and  from  Cj ,  63  £3  will  also  be 

small. 

Now,  considering  z  -  \  and  using  the 
yielding  condition  at  the  plastic -elastic 
front,  we  can  follow  steps  similar  to  those 
that  we  used  for  the  “0”  plane.  As  a  result 
Pii  =  1.11  mm.  It  is  useful  to  find  a  general 
expression  for  e  in  the  “0”  plane.  After 
substitution  of  the  numerical  values  for  a 
and  p,  eqn.  (B12)  becomes,  for  z  =  0, 


1.15 -r 


Aa  AT 


O.SlSr^  +  0.407  -0.885r2 
eg  = - „ - Aa  AT 


_  0.667r3  -  0.575r2  -  0.265 

0.65r2 

e  =  9.077  X10“3  X 


Ao  AT 


/  ,  0.531 

X  (^4.668r2  ~  I0.35r  +  5.952 - - - 

0.42\^2 


This  effective  strain  is  plotted  as  a  function  of 
r  in  Fig.  5.  When  e  =  Cq  =  2  X  10'“*  is  substi¬ 
tuted  in  the  expression  above,  the  theoretical 
plastic  zone  radius  can  be  obtained  from  r  — 
p  =  1.15  mm.  The  radius  of  the  fiber  is  0.5 
mm.  Thus  p  =  1.15^0.5  =  2.3  X  fiber  radius. 
This  radius  includes  the  fiber  size.  Thus  the 
extent  of  the  plastic  zone  measured  from  the 
interface  is  (2.3  —  1)  X  fiber  radius  =  1.3  X 
fiber  radius. 
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Summary 

When  a  metal  matrix  composite  Is  cooled  down  to  room  temperature  from 
the  fabrication  or  annealing  temperature,  residual  stresses  are  Induced  In 
the  composite  due  to  the  mismatch  of  the  thermal  expansion  coefficients 
between  the  matrix  and  whisker. 

An  Investigation  was  undertaken  of  the  extent  of  the  thermal  residual 
stresses  by  an  X-ray  diffraction  technique  as  well  as  of  the  difference  of 
the  yield  stresses  Ao^  between  tension  and  compression  resulting  from  the 

thermal  residual  stresses.  A  theoretical  model  based  on  the  Eshelby's 
method  was  then  constructed  for  the  prediction  of  the  thermal  residual 
str.2sses  and  Aa^.  The  agreement  obtained  was  very  good  between  the 

experimental  results  and  the  theoretical  predictions. 


*  This  research  was  supported  In  part  under  a  contract  from  the  Office  of 
Naval  Research,  Contract  No.  N00014-85K-0007. 
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Introduction 


Metal  matrix  composites  (MMCs) ,  Including  eutectic  composites,  are 
becoming  Important  In  their  application  to  structural  components  which  are 
to  be  used  at  Intermediate  and  high  temperatures.  When  MMCs  are  fabricated 
at  high  temperature  or  annealed  at  a  certain  high  temperature  and  cooled 
down  to  room  temperature,  undesirable  results,  such  as  low  tensile  yield 
stress  and  strength  of  the  MMC  upon  mechanical  testing  at  room  temperature, 
are  produced.  These  results  are  mainly  due  to  residual  stresses  that  are 
caused  by  the  mismatch  of  the  thermal  expansion  coefficients  between  the 
matrix  and  fiber.  This  subject  has  been  studied  by  a  number  of  researchers 
(for  examples  see  references  1-13).  The  residual  stresses  so  Induced  have 

been  observed  In  tungsten  fiber/copper  composltes^^^  and  In  SIC  whlsker/6061 

AT  Most  of  the  models  purposed  to  estimate  the  residual 


AI  composites 


( 1-3  j 

stress  were  based  on  one-dlmenslonal  analysis  ,  continuous  fiber 


system^^’^\  or  spherical  particle  system^ 


The  model  based  on  Eshelby's  equivalent  Inclusion  method^ has  been 


used  to  solve  the  problem  of  thermal  residual  stress'  >  The  advantages 

of  the  Eshelby's  model  are  that  It  can  solve  a  three-dimensional  composite 
system  such  as  short  whisker  composite  and  also  can  take  Into  account  the 


effect  of  the  volume  fraction  of  whisker  (V^)  easily. 


Eshelby's  method  has  been  used  also  for  predicting  the  yield  stress  and 

(15-19) 

work-hardening  rate  of  metal  matrix  composites'  .  The  effect  of  the 

thermally  Induced  residual  stress  on  the  yield  stresses  has  been  discussed 


by  Wakashlma  et  al.'  '  who  predicted  that  the  yield  stress  in  compression 
c  t 

(o^)  exceeds  that  in  tension  (o^)  for  continuous  B  fiber/Al  composites  which 

were  cooled  down,  although  no  comparison  with  the  experimental  data  was 
made. 


In  this  paper  we  focus  on  the  residual  stresses  induced  in  a  short 
whisker  MMC  due  to  the  temperature  drop  and  Its  effect  on  the  yield 
stresses.  The  target  short  whisker  MMC  is  SIC  whlsker/6061  Al  composite. 


Experimental  Procedure 


Materials 


The  target  short  fiber  MMC  was  SIC  whlsker/6061  Al  composite  and  was 
purchased  from  ARCO/SILAG.  The  composite  was  In  the  form  of  an  extruded  rod 


15.5  mm  In  diameter.  Three  different  volume  fractions  of  whisker  (V^)  were 


used:  ■  0,  0.05  and  0.2.  The  composite  was  supplied  as-fabricated  (no 


heat  treatment  made),  was  machined  Into  samples  (Fig.  la  and  lb),  annealed 
for  12  hours  at  810  K,  and  then  furnace  cooled. 


Testing  Methods 


Tension  and  Compression  Tests.  The  testing  was  performed  In  an  Instrora 
testing  machine  using  a  liquid  metal  container  as  the  lower  gripping  device. 
This  method  of  gripping  was  employed  to  ensure  very  good  alignment  of  the 

(20) 


sample.  This  test  procedure  Is  described  In  greater  detail  elsewhere 
Samples  of  two  different  gauge  lengths  were  used  to  determine  If  there  was  a 
gauge  length  effect  (there  was  none).  The  effective  gauge  length  of  the 


Arsenault 


samples  was  determined  by  glueing  strain  gauges  on  the  center  portion  of  the 
sample  and  comparing  the  results  obtained  from  a  clip  on  an  extensometer 
which  mounted  Into  the  "V"  grooves.  If  the  extensometer  is  mounted  directly 
to  the  uniform  gauge  section,  there  is  a  high  probability  that  the  sample 
will  fracture  where  the  '*knife"  edge  of  the  extensometer  makes  contact  with 
the  sample.  Several  tests  in  the  low  stress  range  were  conducted  using  both 
the  extensometer  and  the  strain  gauge;  from  these  tests  the  effective  gauge 
length  was  determined.  Subsequent  tests  were  conducted  using  only  the 
extensomenter.  The  samples,  which  were  tested  in  the  range  from 
-5  -1  -3  -1 

8  X  10  sec  to  2  X  10  sec  ,  showed  no  effect  on  strain  rate. 


Figure  1  -  A  schematic  diagram  of  the 
sample  configuration  which  was  used  in 
the  tension  and  compression  testing. 

The  dimensions  are  In  millimeters. 

Residual  Stress  Measurement  by  X-ray  Diffraction.  The  X-ray 
diffraction  technique  was  used  to  measure  residual  stress.  Two  components 
of  the  residual  stress  were  focused  on:  the  stress  along  the  extrusion 
direction  (longitudinal  direction),  which  coincides  with  the  whisker  axis  of 
a  majority  of  whiskers  (o  ),  and  the  one  in  the  transverse  direction  (o_). 

The  data  were  obtained  and  reduced  by  standard  techniques  as  outlined  by 

Cohen  et  al. 


Experimental  Results 


The  stress-strain  curves  of  .  0,  0.05  and  0.2  composites  are  plotted 

in  Fig.  2a,  2b,  and  2c,  respectively.  The  solid  and  dashed  curves  denote 
the  tensile  and  compressive  test  results,  respectively.  The  yield  stress 
was  measured  as  0.2%  off  set  strain  and  is  indicated  in  the  figures  by 


arrows . 
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Figure  2a  -  The  absolute  stress  vs.  strain  curves  for  tension 
and  compression  test  of  0  volume  %  whisker  material  in  the 

annealed  condition. 
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Figure  2b  -  The  absolute  value  of  stress  vs.  strain  for 
5  volume  %  whisker  composite  in  the  annealed  condition. 
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Figure  2c  -  The  absolute  value  of  the  stress  as  a  function 
of  strain  for  20  volume  %  whisker  material  in  the 
annealed  condition. 

The  residual  stresses  were  measured  by  X-ray  diffraction  on  three  dif¬ 
ferent  values  of  (0,  0.05,  and  0.2)  and  are  tabulated  in  Table  1  where 

the  range  of  scattered  data  is  also  shown.  It  is  noted  that  the  values  of 
and  represent  the  volume  average  quantity  in  the  matrix  (6061  Al).  It 

should  also  be  noted  that  a  compressive  residual  stress  was  found  for  all 
cases. 

Table  I 

Residual  Stresses  X-Ray  Determination 


Sample  Designation 

Trans . 

Long .  o 

MPa 

MPa  ^ 

0 

-37  to  -44 

-18  to  -22 

5 

-17  to  -41.4 

-29  to  -41.4 

20 

-49  to  -57.3 

-24  to  -44.2 

compressive  residual  stress 

Theoretical  Procedure 

The  theoretical  model  used  here  is  based  on  Eshelby's  equivalent 

inclusion  model.  Mori  and  his  co-workers^ extended  Eshelby's  method 
to  predict  the  yield  stress  and  work-hardening  rate  of  aligned  short 

whisker  composites.  Wakashlma  et  al.^^®^  extended  the  above  approach  to 
predict  Oy  by  considering  the  mismatch  of  the  thermal  expansion  coefficients 

of  the  matrix  and  fiber.  Following  the  above  models,  Takao  and  Taya^^^^ 
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have  recently  computed  the  stress  field  in  and  around  a  short  fiber  in  a 
short  fiber  composite  where  the  fiber  is  anisotropic  both  in  stiffness  and 
thermal  expansion. 

In  this  paper  we  focus  on  the  average  thermal  residual  stress  Induced 
In  the  matrix  by  the  cool-down  process  and  also  the  yield  stresses  in  ten- 
t  c 

sion  (Oy)  and  compression  (Oy)  when  the  composite  is  tested  at  the  room 
temperature.  The  former  case  is  essentially  based  on  the  model  by  Takao  and 
Taya^^^\  and  as  the  latter  case,  we  extend  the  model  by  Wakashima  et 
ai.Mo/  account  for  the  bi-linear  stress-strain  curve  of  the  matrix. 

Formulation 

Consider  an  infinite  body  (D)  which  contains  ellipsoidal  whiskers  (11) 
aligned  along  the  x^-axis  (Fig.  3).  This  composite  body  D  is  subjected  to 

the  applied  stress  field  o?..  The  stiffness  tensors  of  the  matrix  (D-fi)  and 

w 

whisker  (fl)  are  denoted  by  C.  and  C.  ,  respectively.  Following 

^  n  1  ( 22) 

Eshelby,  the  transformation  strain^  ^  or  elgenstrain  is  given  in 

the  fiber  domain  II  as 
*  * 
a..,  where  a.,  is  the 

0  ^  J  ^  J 

a-  ■  strain  due  to  the  mis- 

match  of  the  thermal  ex- 

_ T _ T  pansion  coefficients  and 

I  I  the  uniform  plastic 

_ _ _ -  strain  e^j  Is  pre- 

„  scribed  in  the  mat- 

o-a  (C;jKj  .e",)  ,  (15)  .  , 

'  'J  rlx  .  As  far  as  the 

1n\.  stress  field  is  con- 

^  '  w  *  cerned,  the  model  of  Fig. 

^ - ’®ij  ^  ^  equivalent  to  that 

M  //  I  3.  Thus,  the 

VA  y/j  present  problem  is  re- 

y/  ^  duced  to  "inhomogeneous 

yy  Inclusions  problem"  (Fig. 

^  The  model  of 

Fig.  4  will  be  used  not 

§*3  only  to  predict  the  yield 

Xg  stresses  and  work-harden- 

ing  rate,  but  also  to 
*1  compute  the  thermal  resi- 

W  dual  stress.  For  in  the 

I  latter  case,  we  will  set 

/  J  o  p 


.efi) 


i  i 


Figure  3  -  Theoretical  model;  actual  case 
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Figure  A  -  The  equivalent  inclusion  model 
converted  from  Figure  3. 


Following  the  Eshelby's  equivalent  Inclusion  method  modified  for  a 

(27  231 

finite  volume  fraction  of  whiskers'  ’  the  total  stress  field  in  the 
o 


fibers  is  given  by 


Oij  + 


c;..  Je°.  +  e,_.  +  e 


ljk£^®kA 


■k£ 


k£  ■  ^“k£  ■  ®k£^^ 

^ijk£^®k£  ®k£  ®kJl  "  ^®k£  "  ®k£^  "  ®k£^ 


(1) 


where 


o  _  o 

°lj  “  ijk£®k£ 


in  D 


(2) 


•  *.*  ' 
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in  D-R 


(3) 


I1 

• 

1 

K 

t: 

I 


#r,  K , 

•  i 


I 


<a,,> 
ij  m 


In  Eq.  (1),  e,  is  the  fictitious  eigenstrain' 
ke 


which  was  introduced  to 


connect  the  present  problem  to  "inclusion  problem"  and  is  the  average 

strain  disturbance  in  the  matrix  and  related  to  the  average  stress  disturb¬ 
ance  in  the  martix  <o^  ,>  by  Eq .  (3).  <o..>  is  defined  by 

ij  m  ij  m 


<0.  ■> 
ij  m 


^  I 


a.  .dv 


where  ^  is  the  total  volume  of  the  matrix  and  Eqs.  (1) 

and  (4)  is  the  stress  and  strain  disturbance,  respectively,  by  a  single 
fiber  when  it  is  embedded  in  an  infinite  matrix.  Since  the  stress 
disturbance  when  it  is  integrated  in  the  entire  domain  T)  vanishes, 

/  a  dv  =  0  (5) 

D 

The  stress  disturbance  is  obtained  from  Eqs.  (1)  and  (2) 

°ij  ^ijk£^®k£  ®k£  "  ®k£^ 

where  xu.  ^  ^  . 


From  Eqs.  (3),  (5)  and  (6),  we  obtain 


e,  =  -  V  (e.  .  -  e,  .)  (8) 

ij  w"  ij  ij 

■  According  to  Eshelby,  e^^^  is  related  to  the  total  eigenstrain  as 

e .  .  =  S  .  .e,  .  (9) 

ij  ijkil  k£ 

w 

where  ^  is  the  Eshelby's  tensor  and  a  function  of  C  ,  c.  and  the 
IJke  ijk£  ljk£ 

K  geometry  of  the  ellipsoidal  whisker^  .  The  fiber  is  assumed  as  a  prolate 

It  spheroid,  hence  the  whisker  aspect  ratio  t/d  is  a  single  geometrical  para¬ 

meter.  For  simplicity,  we  assume  that  the  matrix  and  fiber  are  isotropic 

w 

I'.  both  in  stiffness  and  thermal  expansion  coefficient.  Thus, 

' are  given  by 

^ljk£  =  ^«lj\£  ^  ^('lk«j£  ^  «i£^j> 

^ljk£  °  ^  ^lj^k£  ^  ^'^lk^j£  '^i£'^k1^ 

“ij  “  -  ^“w  - 

w  w 

In  the  above  equations,  6^^  is  the  Kronecker's  delta,  X  (X  )  and  p  (p  )  are 

Lame's  constants  of  the  matrix  (whisker),  a  (a  )  is  the  thermal  expansion 
I  ” 

coefficient  of  the  matrix  (fiber)  and  AT  is  the  change  in  temperature 
(AT  >  0  corresponds  to  the  temperature  drop). 


The  stress  field  in  the  fiber  is  obtained  from  Eqs.  (6),  (8),  (9)  and 


mm 
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**  ** 


**  ** 


o,  ,  -  (1  -  V  ){(S,  .  e  -  e,  ,  )  X6. ,  +  2u(s,  ,  e  -e,  j} 
ij  w  ‘  kkmn  mn  Tck  ij  Ijmn  mn  ij  ^ 


After  solving  for  in  Eq.  (1)  by  use  of  Eqs.  (8)  and  (9),  we  can  compute 
the  stress  disturbance  in  the  whisker  from  Eq.  (13). 
t  c 

Yield  Stresses  o„,  o„  and  Work-Hardening  Rate 


t  c 

The  method  of  computing  the  yield  stresses  and  and  the  work- 

hardening  rate  of  the  composite  is  described  briefly.  The  total  potential 
energy  of  Fig.  3,  U,  is  given  by 

”  ■  I  ♦  "ij  *  “ij  -  'u  - 

-  /  0°  n  (u°  +  U  +  u  )dS  (14) 

|D|  J  i  1  i 

where  u°,  u^  and  u^  are  the  displacement  components  corresponding  to 
e°  e .  .  and  respectively;  the  index  j  preceded  by  a  comma  denotes  a 

partial  differentiation  with  respect  to  x^;  |D|  is  the  boundary  of  D;  and  n^ 
is  the  j-th  component  of  an  unit  vector  outer  normal  to  |d|.  Then,  the 
change  in  u  due  to  the  change  in  the  plastic  strain  6  e^^  is  given  by 

6U  =  -  /  (o°j  +  Oj^j)6  eJjdV  (15) 

In  the  above  deviation,  the  Gauss'  divergence  theory  and  the  following 
equation  were  used 

+  SUjj)dV  .  0 

Noting  that  the  plastic  strain  exists  only  in  the  matrix  (D-O),  Eq.  (15)  is 
reduced  to 

6U  =  -  dej.  {(1-V  )a°  -Vo,,}  (16) 

ij  ^  w'  ij  w  ijJ 


Under  the 


uniaxial  stress  along  the  x^-axls  (o  ),  a°  ®nd  ®re  given  by 

^  o  ij  ij 


-l/2e 

1 

-l/2e 


where  the  six  components  of  and  ejj  are  expressed  in  the  order  of  (ij)  ■ 

11,  22,  33,  23,  31,  and  12,  and  Op  is  the  plastic  strain  along  the  x^- 

axis.  On  the  other  hand,  the  energy  dissipation  due  to  the  plastic  work  in 
the  matrix,  6Q,  is  given  by 

6Q  -  (1-V  )a  6e  (18) 

w  y  p 


. '  .*•  V-  -  •  . 
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where  a  is  the  flow  stress  of  the  matrix  for  the  bilinear  model 

y 

Oy  -  0°  +  E^(e  -  Cp)  (19) 

where  0°  are  the  initial  yield  stress  and  tangent,  modulus  of  the 

matrix,  respectively,  and  e  is  the  total  strain.  Since  6U  +  60  »  0,  we 
obtain 

"o  ■  °y  (PV;;)'  <”33  - 

In  the  above  deviation,  ^22^  e^^  =  6^2  “  “  ®p  used.  Combining 

the  solutions  of  o  in  Eq.  (13)  and  Eq.  (20),  we  can  obtain  the  yield 
^  J  t 

stress  of  the  composite  in  tension  (o^  =  o^)  and  that  in  compression 
(Oy  -  o^)  as 

=  ^0  ^1  “m^"^ 

c  o  ^21) 

a  /a  •  -  C,  a  AT 

y  y  0  1  m 

Similarly,  the  work-hardening  rate  of  the  composite,  E^  is  obtained  as 

ll/E,  -  Cj  *  C3  (i)  (22) 

In  the  above  equations,  Cq,  Cp  C2 ,  and  are  functions  of  the  mechanical 

properties  of  the  matrix  and  whisker  and  the  fiber  aspect  i/d,  E  is  the 
Young's  modulus  of  the  matrix. 

Thermal  Residual  Stress 

When  the  composite  is  cooled  down  by  AT,  the  thermal  residual  stress  is 
Induced  in  the  composite.  The  theoretical  model  for  this  problem  is  the 

same  as  that  shown  in  Fig.  4  except  that  =  e^j  »  0  in  the  present 

case.  Hence,  the  formulation  up  to  Eq.  (13)  is  valid  and  will  be  used  to 
compute  the  thermal  residual  stress.  Once  the  stress  within  the  whisker 

is  computed,  the  stresses  just  outside  the  whisker  can  be 

obtained  by  the  following  relation^ 


(out) 


^  ^pqmn^  Stilj®lj”i”i 


IjVn  v(X+2w) 


+  (23) 

mn 

(in)  ** 

where  is  given  by  Eq .  (13),  e^^^  is  solved  by  Eq.  (1),  \  and  p  are 

previously  defined  and  n^  is  the  1-th  component  of  the  unit  vector  pre- 

pendlcular  to  the  surface  of  the  whisker.  If  one  wants  to  compute  the 
stresses  just  outside  the  equator  of  the  ellipsoidal  whisker  where  n  -  (1, 

0,  0),  we  can  obtain  the  stresses  in  polar  coordinates,  o  ,  o  and  o  there 

re  z 

as 


.  ^  *- 


-  > 

:  •/  • 
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(out)  -(in) 

Or  -  Or 


(out)  (In)  -  f  1  **  V  **, 

"e  -“e  *  Ill'll -"T^Ssl 


(24) 


(out)  (In)  «  r  V  **  1  **, 

-‘’z  +3^633} 

where  v  is  the  Poisson's  ratio  of  the  matrix. 


Next  the  average  stress  field  in  the  matrix,  <0  >  ,  can  be  also 

obtained  from  Eqs.  (3),  (8)  and  (9)  as  ^  ” 


<0  > 

Ij  m 


-  V  C,,,  .(S,  ,  .. 

w  IjkA  kimn  mn 


**  ** 

e  -  e,..) 


'ki' 


(25) 


The  average  stress  In  the  matrix,  <0  >  ,  <o.>  and  <0  >  ,  will  be  computed 

r  m  6  m  z  m 


by  setting  Ij  =  11,  22,  and  33  in  Eq.  (25).  It  should  be  noted  here  that 


<0  > 
r  m 


<a„>  due  to  the  assumptions  of  an  aligned  short  whisker  composite 
6  m 


which  results  in  the  transverse  Isotropy  of  the  volume  average  quantity. 

Theoretical  Results 


The  thermo-mechanical  data  of  the  matrix  and  whisker  for  the  theore¬ 
tical  calculations  are  obtained  from  the  stress-strain  curve  of  the  matrix 
(Fig.  2a)  and  the  material  properties  handbook. 


Annealed  6061  A1  matrix  (bilinear  model); 
E  -  47.5  GPa 


=  2.3  GPa 


V 

a 


47.5  MPa 
10" 


(26) 


0.33 
23.6  X 


V°C 


SIC  Whisker! 


427  GPa 
0.17 


a^_  =  4.3  X  10"^/°C 


(27) 


ijd  =  1.8 

where  the  average  value  of  the  fiber  aspect  ratio  (l/d)  was  used 
bilinear  stress  strain  curve  of  the  matrix  Is  Indicated  by  a  dash-dot  line 
In  Fig.  2a.  The  temperature  drop  AT  Is  defined  as 


and  the 


AT 


T  -  T 

h  0 


(28) 


where  T^  Is  taken  as  the  temperature  below  which  dislocation  generation  is 


minimal  during  the  cooling  process and  Tq  is  the  room  temperature, 
for  the  present  composite  system  AT  is  set  equal  to  200  K. 


Thus , 


c  ( 

The  theoretically  predicted  o  and  0  can  be  obtained  from  Eq.  (21). 

t^  c  ^t 

o  (o  -  o  )  for  various  V„  are  shown  as  the 

y  y  y  t 

solid  line  in  Fig.  5.  The  differences  between  o  and  o  increase  with 

y  y  c  t 

increasing  V  .  The  experimentally  determined  differences  between  o  and  o 
w  y  y 


The  differences  between  a  and 

y 


are  represented  by  the  open  circles  in  Fig.  5.  The  experimental  values  of 
the  yield  stresses  are  taken  as  the  stresses  0.2%  off-set  strain.  It 
follows  (Fig.  5)  that  good  agreement  is  obtained  between  the  experimental 


ii  r  .  'v  “  O  'y  V  •  •  ' 
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Figure  5  -  The  difference  in  yield  stress  between 
compression  and  tensin  as  a  function  of  volume  % 
silicon  carbide  whisker. 


and  theoretical  results  of  the  difference  in  the  compressive  tensile  yield 

c 

stresses.  However,  the  experimental  values  of  a  are  greater  than  the 

c  ^  t 

theoretical  values  of  o  ,  and  similarly  the  experimental  values  of  o  are 

^  t  ^ 

greater  than  the  theoretical  values  of  o  . 

y 

From  the  data  given  by  Eqs.  (26)  -  (27)  and  the  use  of  Eq.  (25),  we 
have  computed  the  average  stresses  in  the  matrix  and  the  stresses  just 
outside  the  fiber  and  plotted  schematically  the  and  along  the  x  (or 

X2)  and  Xj  axes  in  Fig.  6. 
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Figure  6  -  Is  a  schematic  of  the  stress  distribution 
in  the  matrix  and  in  the  reinforcement  due  to  the 
difference  in  thermal  coefficient  of  expansion 
between  the  silicon-carbide  and  aluminum. 


Next,  the  thermal  residual  stresses  averaged  in  the  matrix  of  the  SiC 

whi8ker/6061  A1  are  predicted  by  Eq.  (25)  and  the  results  on  <o  >  and 

T  m 

<0, >  are  plotted  in  Fig.  7  as  a  function  of  the  volume  fraction  of  whisker 
Li  m 

(V^),  where  the  subscripts,  T  and  L  denote  the  component  along  the  trans¬ 
verse  direction  (r  and  6)  and  longitudinal  direction  (z).  The  average 
theoretical  thermal-residual  stress  is  predicted  to  be  tensile  in  nature. 


Arsenault 


Figure  7  -  Predicted  residual  stress  In  the 
matrix  for  the  transverse  and  longitudinal  directions. 

Discussion  and  Conclusions 

The  theoretical  predictions  and  experimental  results  are  in  very  good 
agreement  In  most  cases,  and  In  some  cases  the  differences  are  to  be 
expected. 

The  theoretical  model  Is  based  on  an  extension  of  previous  work  by 

Eshelby^^^\  Mura  and  Taya^^^\  and  Tanaka  and  Morl^^^\  and  it  can  predict 
the  yield  stress  in  tension  and  compression  and  the  thermal  residual  stress. 
The  predicted  values  of  the  yield  stress  in  tension  and  compression  are  less 
than  the  experimentally  determined  values  of  the  yield  stress  in  tensional 
compression.  This  difference  between  the  theoretical  and  experimental 
stresses  is  due  to  the  fact  that  the  V„  ■  0  curves  of  stress  vs.  strain  were 

used  for  the  matrix  in  the  composite  cases.  It  has  been  demonstrated  by 

Arsenault  and  Fisher^®^  and  Vogelsang,  et  al.^^^  that  there  is  a  much  higher 
dislocation  density  in  the  annealed  composite  matrix  than  in  the  ~  0 

material.  Therefore,  the  matrix  is  stronger  then  the  annealed  ■  0  matrix 

alloy.  However,  this  Increase  in  matrix  strength  due  to  a  higher  dlsloca- 


i 


V. 
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tlon  density  does  not  influence  the  difference  in  yield  stress  in  the  com¬ 
pression  vs.  tension.  The  theoretical  prediction  is  that  the  yield  stress 
in  compression  should  be  higher  than  the  yield  stress  in  tension,  and  this 
is  exactly  what  is  observed  experimentally. 

The  theoretical  prediction  that  the  compressive  yield  stress  is  higher 
than  the  tensile  yield  stress  suggests  that  a  tensile  residual  stress  exists 
in  the  matrix  of  composite.  An  average  tensile  residual  stress  is  predicted 
in  the  i^atrix  as  shown  in  Figs.  6  and  7.  However,  the  experimentally  deter¬ 
mined  resldi'dl  stress  is  compressive,  l.e.,  opposite  of  the  prediction. 

This  difference  can  qualitatively  be  explained  in  terms  of  the  magnitude  of 
the  tensile  and  compressive  residual  stresses  in  the  matrix.  If  the  com¬ 
pressive  regions  of  the  matrix  were  eliminated  and  only  the  tensile  regions 
were  examined,  then  there  would  be  tw  change  in  the  diffraction  peak  condi¬ 
tions,  because  the  limit  of  detection  is  50  MPa.  If  a  similar  experiment 
were  conducted,  but  in  this  case  only  the  compressive  regions  were  kept  then 
there  would  be  a  change  in  peak  position  because  the  compressive  residual 
stress  is  greater  than  the  detection  limit.  Now,  if  both  tensile  and  com¬ 
pressive  regions  are  combined,  the  net  result,  a  small  compressive  residual, 
will  be  experimentally  detected,  which  is  in  agreement  with  the  X-ray 
residual  stress  measurements. 

It  is  possible  to  arrive  at  the  following  conclusions  from  the 
theoretical  model  and  the  experimental  results. 

1.  The  theoretical  model  predicts  a  higher  compression  yield  stress  than 
tensile  yield  stress,  which  is  in  agreement  with  the  experimental 
results . 

2.  The  absolute  magnitude  of  the  predicted  yield  stress  (tensile  and 
compressive)  is  less  than  the  experimental  yield  stress  due  to  an 
Increase  in  matrix  strengthening.  The  increased  matrix  strength  is  due 
to  a  higher  dislocation  density  in  the  matrix;  the  higher  dislocation 
density  is  a  consequence  of  the  difference  in  coefficient  of  thermal 
expansion  between  the  A1  alloy  matrix  and  the  SIC. 

3.  The  magnitude  of  the  average  residual  stress  is  small,  but  at  the  SiC-Al 
alloy  interface  there  can  be  a  large  compressive  stress. 

4.  The  residual  stress,  in  the  region  between  SIC  whiskers,  is  tensile  and 
this  is  the  likely  region  where  plastic  deformation  would  begin.  For  in 
this  region  the  matrix  contains  a  lower  dislocation  density  than 
adjacent  to  the  SIC  whisker,  i.e.,  in  this  region  the  matrix  is  weaker. 
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Abstract 


An  investigation  was  undertaken  to  evaluate  the  strength  o"^ 
the  interfacial  bond  between  SiC  and  the  6061  A1  alloy  matrix. 
Approximate  bounding  analysis  provided  the  conditions  for  separa¬ 
tion  of  the  inclusion  from  the  ductile  matrix  under  the  local 
negative  pressure  (triaxial  tensile  stress).  The  experimental 
data  was  analyzed  to  determine  the  local  interfacial  stresses  at 
the  particle-matrix  interfaces.  The  lower  bound  value  of  the 
bona  strength  was  determined  to  be  equal  to  at  least  1690  I'lpa. 
Debonding  of  SiC  particulates  from  the  matrix  was  found  to  be  a 
rare  event. 


Introduction 


The  quality  of  interfacial  bond  between  A1  and  SiC  is  impor¬ 
tant  in  composite  strengthening,  independent  of  the  strengthening 
mechanism  that  is  assumed  to  be  operative.  One  of  the  key 
elements  in  the  continuum  mechanics  treatment  of  the  composite 
strengthening  is  that  the  interface  transfers  the  load  from  the 
matrix  into  the  reinforcement  1- ^  ^ .  The  rule  of  mixtures  that 

emerges  from  the  continuum  mechanics  relates  a  given  property  o"*^ 
the  composite  to  the  properties  of  its  constituent  materials. 

The  credibility  of  this  approach  depends  on,  among  other  things, 
the  bond  integrity  and  the  efficiency  of  the  load  transfer  from 


the  matrix  to  the  reinforcement. 


2 


Recently  it  has  been  shown!- that  in  Al/ SiC  systems  the 
interfaces  are  a  major  source  of  the  dislocation  generation  (on 
cooling  from  fabricating  temperature  due  to  the  difference  in  the 
coefficient  of  thermal  expansion).  The  resultant  dislocation 
densities  at  the  Al/SiC  interfaces  can  be  very  high  (10^^  to  10^'^ 
cm~^),  which  significantly  contributes  to  the  overall  composite 
strength.  Interestingly,  when  an  Al/SiC  system  is  subject  to 
heating,  the  Al/SiC  interface  also  generates  dislocations  (pro¬ 
viding  that  Al/SiC  bond  is  quite  st rong)  ^ .  If  the  bond  is 

weak,  the  Al  just  pulls  away  from  the  particle  and  is  free  to 
expand  without  any  restrictions  on  the  part  of  the  SiC,  and  in 
the  subsequent  cooling  no  dislocations  are  generated.  Therefore, 
a  good  bond  is  required  to  produce  the  high  dislocation  density. 

Several  investigators  have  reported  that  the  bond  between 
Al/SiC  is  generally  good!-®’9j,  but  a  systematic  evaluation  o^  the 
bond  strength  has  not  been,  to  our  knowledge,  reported  in  the 
literature.  Direct  measurements  of  the  interfacial  shear 
strength  in  metal  matrix  systems  have  been  made  using  flat  plate 
and  fiber  pullout  tests!- The  fiber  pullout  test,  which  is  cf" 
more  interest  for  this  investigation,  was  used  to  determine 
interfacial  shear  strength  in  Cu/W  and  Cu/I'^o  systems!- ^  ^  . 

Several  attempts  were  made  to  carry  out  a  direct  measurement  o-f 
the  interfacial  bond  shear  strength  in  Al/SiC  using  the  pullout 
test  Que  to  its  apparent  simplicity.  However,  all  these  attemnts 
resuxtea  in  brittle  failure  of  SiC  (single  crystal,  very  large 
grain,  and  sintered  SiC  were  used)  at  very  moderate  loads. 


3 

In  addition,  an  effort  was  made  to  perform  a  punch  test 
using  Ai  disks  with  pieces  of  the  SiC  embedded  in  the  center. 

This  effort  also  resulted  in  premature  failure  of  the  SiC. 

Another  way  to  evaluate  the  interfacial  bond  strength  is 
based  on  the  general  understanding  of  the  mechanism  of  the 
ductile  f racture  1- ^  ^  5  ]  applied  for  the  determination  o"^ 

the  interfacial  bond  strength  in  spheroidized  104^  steel,  Cu-0.6 
pet.  Cr  alloy,  and  managing  steel  containing  Fe-^C.  Cu-Cr  and  '^iC 
particles,  respectively  I- ^  ^ ^  •  Providing  that  th'^  ductile 
fracture  starts  by  void  nucleation  at  the  inclusions  and  second 
phase  particles,  the  theoretical  bounding  analysis-''-!  allows  the 
evaluation  of  the  interfacial  strength  in  the  systems  with  the 
small  volume  fraction  of  the  second  phase. 

The  purpose  of  this  investigation  was  to  experimentally 
determine  the  bond  strength  between  SiC  and  Al  in  a  commercial 
composite  material  using  the  analysis  of  separation  o-^  the 
inclusion  from  the  matrix  under  the  negative  pressure. 


Experimental  Procedure 


One  volume  percent  (IV'^)  SiC  particulate  in  a  6061  Al  alloy 
matrix  composite  purchased  from  DWA  was  used  for  this  investiga¬ 
tion.  The  low  volume  percent  of  the  particulate  was  necessary  in 
order  to:  1)  provide  favorable  conditions  for  a  ductile  fracture 
and  thus  activate  the  mechanism  of  void  nucleation  and  growth 
type  failure,  2)  satisfy  criteria  of  non-interacting 
particles!- ^  ^ .  At  higher  volume  concentrations,  the  Al/SiC 
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composites  exhibited  very  little  ductility.  The  theoretical 
analysis  of  the  inclusion  separation  from  the  matrix  that  was 
used  in  this  investigation  was  based  on  the  assumption  that  there 
is  no  interaction  between  the  particles,  which  is  a  reasonable 
statement  for  a  small-volume  concentration. 

The  composite  material  was  machined  to  a  standard  tensile 

test  specimens.  In  order  to  introduce  a  local  triaxial  stress 

state,  a  circumferential  groove  was  electrical  discharge  machined 

(EDM)  in  the  center  of  the  specimen  (Eig.  1).  Subsequently,  the 

specimen  were  solution-annealed  at  823  K  for  12  hours  and  jce 

water  quenched  in  order  to  keep  the  Mg2Si  phase  in  the  solution 

and  thus  limit  the  presence  of  particles  other  than  SiC.  After 

quenching,  the  samples  were  placed  into  the  freezer  of  a 

commercial  refrigerator  where  they  were  stored  at  265  K  and 

within  a  few  minutes  the  samples  were  able  to  warm  up  to  the 

testing  temperature  of  295  K.  The  specimens  were  tested  in 

tension  to  fracture  using  the  Instron  testing  machine  at  cross- 

cm 

head  velocity  of  0.05  ^Yn'  fractured  halves  were  EDM,  cut 

longitudinallp  (parallel  to  the  tensile  axis),  and  mechanically 
polished.  Final  polishing  was  done  using  a  colloidal  silica 
suspension  as  abrasive  (Buehler's  Mastermet)  on  the  special 
moderate  length  nap  cloth  (Buehler's  Mastertex).  Since  this 
suspension  has  a  slightly  alkali  reaction,  the  final  polishing 
was  accomplished  both  mechanically  and  chemically.  Therefore, 
the  surface  was  ready  for  inspection  without  the  necessity  of 
etching  to  delineate  the  SiC  particulates.  Polished  sections 
were  examined  for  voids  and  particles  separation  using  the 


scanning  electron  microscope  (SEM). 


Experimental  Results  and  Discussion 

The  experimental  technique  adopted  in  this  work  was  based  on 

the  determination  of  local  stresses  during  plastic  deformation 

under  a  triaxial  stress  state^-^^^.  The  bonding  analysis- 

showed  that  the  interfacial  stress  can  be  expressed  as 

r  r 

=  oj  +  Y(¥P)  (1  ) 

where  Ofj  is  the  local  negative  pressure  (triaxial  tensile 
stress),  Y(eP)  is  true  flow  stress  in  tension  corresponding  to 
the  local  average  plastic  strain,  had  the  second  phase  particle 
been  absent.  This  analysis  was  based  on  the  assumption  that: 

(1)  particles  have  an  equiaxed  shape,  and  (2)  volume  fraction  of 
the  second  phase  is  small.  Also,  it  can  be  mentioned  that  SiC  is 
rigid  and  undeformed.  This  physical  property  agrees  very  well 
with  the  assumption  of  rigidity  of  the  inclusions'^  ^  ^  ^ . 

The  distribution  of  the  triaxial  tensile  stress  along  the 
radial  line  in  the  plane  of  the  groove  was  obtained  from  the 
theory  of  stress  concentrations  on  circumferentially  grooved 
elastic  bars  originally  considered  by  Neuber  and  then 
expressed^- '  ^  ^ 


as : 


6 


0  j 

where  —  represents  the  triaxiality,  om  is  the  negative  pressure, 
°o 

is  the  flow  stress  or  average  ligament  stress,  z  is  equal  to 
the  vertical  distance  along  the  z  axis,  i.e.,  tensile  axis  center 
line  of  sample,  and  r  is  the  distance  from  the  z  axis.  Since 
triaxiality  a^/o^  reaches  its  maximum  value  at  the  bottom  of  the 
groove,  r  should  be  equal  to  a,  where  2a  is  the  diameter  of  the 
ligament  (see  Fig.  1).  Parameters  c  and  a  are  defined  as 
follows : 


1  +  ^  +  /I  +  I 
2(2  +  I  +  ✓I  + 


(?) 


where  a  is  the  radius  of  the  ligament  and  R  is  the  radius  o-f  the 
groove  (see  Fig.  lb).  The  longitudinal  sections  of  the  tested 
specimen  are  shown  on  Fig.  2.  In  general,  the  number  of  voids 
associated  with  debonding  of  SiC  particulates  was  much  smaller 


than  the  total  number  of  voias  related  to  the  fracture.  Several 


examples  of  the  areas  where  debonded  SiC  particulates  can  be 
observed  are  shown  in  Fig.  5.  Debonding  shown  in  Fig.  3  is  a 
rather  rare  event  and  is  not  typical  of  Al/SiCp  fractured 
samplesl^  ^  ®  ^ . 

According  to  the  Neuber  analysis,  the  largest  negative 
pressure  occurs  at  the  outer  surface  of  the  groove  (see  Fig. 
1).  Since,  in  general,  no  separation  between  the  SiC  partic¬ 
ulates  and  A1  matrix  was  observed,  it  is  reasonable  to  assume 
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that  for  Al/ SiC  bond 


„  .  „mSiX 

°rr  ^  °T 


+  Y(eP) 


value  of  Y(e^)  was  obtained  from  stress-strain  curve 


as 


r  ^ 

Y(eP)  =  ^  =  o  =  flow  stress 

A  ^  O 


(5) 


where  will  equal  the  load  at  failure  which  will  equal  750  kg 
and  kf  will  equal  the  area  across  the  grooved  region.  The  values 
of  c  and  a  we  get  from  expressions  (3)  and  (4),  respectively.  At 
the  bottom  of  the  circumferential  grove  (z  =  0)  these  values  are 
c  =  0.48  and  a  =  2.3  mm. 

Substituting  Eqs.  (5)  and  (2)  into  Eq .  (!)  we  obtain: 


O  =0 

rr  o 


/1  -  {^) 


^^2 


+  0  =1690  MPa 

0 


Thus,  the  lower  bound  value  for  Al/SiC  bond  strength  is  16o0  MPj 
the  Neuber  analysis  is  applicable  to  a  linear  solid,  i.e.,  to  a 
material  with  linear  strain  hardening.  Prom  the  stress-strain 
curves  obtained,  it  was  felt  that  the  behavior  of  the  tested 
specimen  could  be  approximated  to  a  linear  body. 


Conclusions 


Providing  that  the  above-mentioned  assumptions  are  reason¬ 
able  (and  there  is  no  reason  to  question  the  assumptions  if  a 
lower  limit  of  the  bond  strength  is  determined),  it  can  be 
concluded  that  the  6061  Al/SiCp  bond  strength  is  at  least  16P0 
MPa.  (This  is  40  times  higher  than  the  yield  stress  of  the 
annealed  6061  A1  alloy.)  It  is  thought  that  debonding  can  be 
attributed  to  the  defects  during  fabrication  of  the  composite 
material.  Ductile  fracture  of  low-volume  percent  6061 /SiCp 
occurs  by  void  nucleation  at  SiC  particulates  and  also  at  the 
matrix  imperfections  (preexisting  voids,  inclusions,  MgSij, 
etc.).  A  separate  investigation,  however,  is  required  to 
determine  explicitly  the  origin  of  the  ductile  fracture  in  Al/SiC 
systems.  It  is  realized  that  the  assumption  of  the  equiaxiality 
of  the  SiCp  is  rather  poor.  This  results  in  further  under¬ 
estimation  of  the  bond  strength  since,  under  the  same  macro 
stress  state,  irregular  shaped  particles  would  have  higher  stress 
concentrations,  and,  thus,  separation  from  the  matrix  would  be 
expected  to  occur  at  the  earlier  stage  compared  to  a  spherical 
particle.  If  the  irregular  shape  of  the  particle  could  be  taken 
into  account,  the  lower  limit  would  increase. 
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List  of  Figures 

Fig.  1  A  schematic  view  of  the  specimen  (a)  and  enlarged 

portion  of  the  groove  (b)  showing  groove  geometry  and 
location  of  the  maximum  triaxial  stresses. 

Fig.  2  Electron  micrographs  of  the  different  areas  of  the 

longitudinal  cross  sections  of  the  fractured  1V^  Al/SiCp 
specimens . 

(a,b)  -  taken  in  the  vicinity  of  the  bottom  of  the 
groove.  No  apparent  debonding  is  present  at  the 
position  of  maximum  triaxial  stress. 

(c,d,e,f,g)  -  taken  below  the  fracture  surface  away  from 
the  groove. 

(h,i,j)  -  taken  far  away  from  the  fracture  and  represent 
the  bulk  of  the  specimen. 

Fig.  3  Electron  micrographs  of  the  areas  below  fracture  surface 
at  the  bottom  of  the  groove  (a)  and  away  from  it  (b-e) 
where  debonding  was  observed. 
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Electron  micrographs  of  the  different  areas  of  the 
longitudinal  cross  sections  of  the  fractured  IVf-  Al/EiCp 
specimens. 

(a,b)  -  taken  in  the  vicinity  o!'  tne  bottom  ol  tne 
groove.  No  apparent  debonding  is  present  at  the 
position  of  maximum  triaxial  stress. 

(c,d,e,f,g)  -  taken  below  the  fracture  surface  away  "'ror, 
the  groove. 


(h,i,j)  -  taken  far  away  from  th*-  ‘‘racture  and  represent 


the  bulk  of  the  specimen 


at  the  bottom  of  the  groove  (a)  and  away  from  it  (b-e) 
where  debonding  was  observed. 
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%:  An  In  Situ  HVEM  Study  of  Dislocation  Generation 
■'  at  Al/SiC  Interfaces  in  Metal  Matrix  Composites 

■  MARY  VOGELSANG,  R.  J.  ARSENAULT,  and  R.  M.  HSHER 


IT 


Annealed  aluminum/silicon  carbide  (Al/SiC)  composites  exhibit  a  relatively  high  density  of  dis¬ 
locations,  which  are  frequently  decorated  with  fine  precipitates,  in  the  A1  matrix.  This  high  dislocation 
density  is  the  major  reason  for  the  unexpected  strength  of  these  composite  materials.  The  large 
difference  (10:1)  between  the  coefficients  of  thermal  expansion  (CTE)  of  A1  and  SiC  results  in 
sufficient  stress  to  generate  dislocations  at  the  Al/SiC  interface  during  cooling.  In  this  in  situ 
investigation,  we  observed  this  dislocation  generation  process  during  cooling  from  annealing  tem¬ 
peratures  using  a  High  Voltage  Electron  Microscope  (HVEM)  equipped  with  a  double  tilt  heating 
stage.  Two  types  of  bulk  annealed  composites  were  examined:  one  with  SiC  of  discontinuous  whisker 
morphology  and  one  of  platelet  morphology.  In  addition,  control  samples  with  zero  volume  percent 
were  examined.  Both  types  of  composites  showed  the  generation  of  dislocations  at  the  Al/SiC 
interface  resulting  in  densities  of  at  least  10'^  m'^.  One  sample  viewed  end-on  to  the  whiskers  showed 
only  a  rearrangement  of  dislocations,  whereas,  the  same  material  when  sectioned  so  that  the  lengths 
of  whiskers  were  in  the  plane  of  the  foil,  showed  the  generation  of  dislocations  at  the  ends  of  the 
whiskers  on  cooling.  The  control  samples  did  not  show  the  generation  of  dislocations  on  cooling 
except  at  a  few  large  precipitate  particles.  The  results  support  the  hypothesis  that  the  high  dislocation 
density  observed  in  aimealed  composite  materials  is  a  result  of  differential  thermal  contraction  of  A1 
and  SiC.  The  SiC  particles  act  as  dislocation  sources  during  cooling  from  annealing  temperatures 
resulting  in  high  dislocation  densities  which  strengthen  the  material. 
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I.  INTRODUCTION 

The  incorporation  of  20  vol  pet  discontinuous  SiC  whisk¬ 
ers  into  a  6061  A1  matrix  increases  the  yield  strength  of 
annealed  powder  compacted  6061  A1  alloy  by  more  than  a 
factor  of  two.  This  increase  in  strength  cannot  be  explained 
directly  by  continuum  mechanics  theories.  Continuum  me¬ 
chanics  formulations  developed  by  Piggott'  and  applied  to 
the  case  of  discontinuous  Al/SiC  composites  by  Arsenault’ 
predict  an  ultimate  strength  of  only  186  MPa  for  20  vol  pet 
SiC  composite,  whereas  the  measured  value  of  ultimate 
strength  for  this  material  is  448  MPa.  Arsenault  and  Fisher’ 
proposed  that  the  increased  strength  could  be  accounted  for 
by  a  high  dislocation  density  in  the  A1  matrix  which  is 
observed  in  bulk  composite  material  annealed  for  as  long  as 
12  hours  at  810  K. 

The  dislocation  generation  mechanism  proposed  by 
Arsenault  and  Fisher  to  account  for  this  high  dislocation 
density  is  based  on  the  large  difference  ( 10 : 1 )  in  coefficients 
of  thermal  expansion  (CTE)  of  A1  and  SiC.‘  When  the 
composite  is  cooled  from  elevated  temperatures  of  an¬ 
nealing  or  processing,  misfit  strains  occur  due  to  differential 
thermal  contraction  at  the  Al/SiC  interface  which  are  suf¬ 
ficient  to  generate  dislocations. 

Chawla  and  Metzger,  in  an  elegant  investigation  of  Cu/W 
composites  using  etch-pitting  techniques,  observed  a  high 
dislocation  density  at  the  Cu/W  interface  which  decreased 
with  increasing  distance  from  the  interface.’  They  observed 
that  if  the  volume  fraction  of  W  was  15  pet,  the  minimum 
dislocation  density  in  the  matrix  was  7  x  10”  m  ’  in- 
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creasing  to  4  x  10”  m*’  at  the  interface  of  W  and  Cu,  and 
concluded  that  the  dislocations  were  caused  by  the  differ¬ 
ences  (4: 1)  in  CTE  of  Cu  and  W.  Recalling  that  the  CTE 
difference  between  A1  and  Si  is  10: 1 ,  i.e.,  more  than  twice 
as  great  as  the  Cu/W  system,  one  would  expect  thermal 
stresses  in  Al/SiC  to  be  sufficient  to  generate  dislocations  in 
this  composite. 

Other  causes  may  also  contribute  to  the  high  dislocation 
density  observed  in  annealed  Al/SiC  material.  Dislocations 
are  introduced  into  this  material  during  the  plastic  defor¬ 
mation  processes  of  manufacturing,  such  as  extrusion.  Dur¬ 
ing  annealing,  the  dislocations  introduced  during  processing 
may  not  be  annihilated;  they  could  be  trapped  by  the  SiC, 
resulting  in  a  high  dislocation  density  after  aiuiealing. 

It  is  important  to  determine  the  origins  of  the  high  dis¬ 
location  density  in  the  composite  since  the  strength  of  the 
composite  depends  on  the  high  density.  If  the  differen¬ 
tial  thermal  contraction  is  the  cause  of  the  dislocations,  as 
Arsenault  and  Fisher’  suggest,  then  dislocations  should  be 
observed  being  generated  in  a  composite  thin  foil  sample 
on  cooling  from  annealing  temperatures  in  an  in  situ, 
HVEM  experiment. 

In  situ  dynamic  HVEM  experiments  have  certain  advan¬ 
tages  over  other  experimental  techniques.  The  major  advan¬ 
tage  is  that  direct  observation  of  a  dynamic  process  altering 
a  microstructure  is  possible  while  Ae  deforming  force,  in 
this  case  a  thermal  stress,  is  operating.  Operating  at  higher 
voltages  allows  penetration  of  thicker  samples  so  that  sur¬ 
face  effects  are  minimized  and  bulk  behavior  is  more  closely 
approximated.  Also,  a  high  voltage  microscope  can  better 
accommodate  special  stages  required  for  in  situ  work  be¬ 
cause  of  the  large  pole  piece  region. 

Numerous  in  situ  HVEM  heating  stage  investigations  of 
A1  have  been  performed.  Hale  etal.*  and  Caillard  and 
Martin’  investigated  dislocation  motion  during  creep  at 


elevated  temperatures  using  HVEM.  Kivilahti  et  al}  ob¬ 
served  an  Al-2  pet  Mg  alloy  in  situ  during  recovery 
processes  at  elevated  temperatures  recording  dislocation 
interactions  on  videotape.  Shimotomori  and  Hasiguti’  ob¬ 
served  in  situ  prismatic  punching  of  dislocation  loops  at 
second  phase  precipitates  in  an  Al-1 .3  pet  Li  alloy.  Electron 
irradiation  of  Al  at  elevated  temperatures  has  b^n  exten¬ 
sively  studied  using  HVEM.'® " 

There  have  been  several  non-dynamic,  non-in  situ  TEM 
investigations  of  dislocations  about  particles  in  a  metal  ma¬ 
trix.  Weatherly'"'  observed  multiple  slip  mode  dislocation 
tangles  around  silica  in  Cu,  and  concluded  they  were  caused 
by  differential  thermal  contraction  of  the  two  materials  on 
quenching.  Ashby  et  alV  observed  dislocations  about  pres¬ 
surized  silica-Cu,  noting  a  critical  size  dependence  for  dis¬ 
location  generation.  Williams  and  Garmong'^  reported  a 
high  incidence  of  dislocations  at  the  Ni/W  interface  in  this 
directionally  solidified  eutectic  composite. 

Calculations  of  the  dislocation  density  in  Al/SiC  due  to 
thermal  stresses  predict  high  dislocation  densities.  The  mis¬ 
fit  strain  which  develops  at  the  circumference  of  a  1  ^m 
diameter  SiC  particle  due  to  differential  thermal  contraction 
during  cooling  is  approximately  1  pet.  The  plastic  strain  at 
one-half  the  interparticle  spacing,  obtained  from  Lee 
et  al.,'^  ranges  from  1  to  2  pet.  The  dislocation  density  can 
be  simply  calculated  from  the  following  equation: 

Ep  =  pLb.  [11 

where  e,  is  the  plastic  strain  (1  pet),  p  is  the  dislocation 
density  (m’’)  generated,  L  is  the  average  distance  moved  by 
the  generated  dislocations,  which  was  taken  to  be  1/2  the 
inter-whisker  spacing,  i.e.,  2  pm,  and  b  is  the  Burgers 
vector  of  Al.  The  p  obtained  is  1.8  x  10'*  m"^ 
Consideration  of  another  type  of  dislocation  described  by 
Ashby'®  predicts  additional  dislocations  in  the  matrix.  These 
dislocations  are  called  ''geometrically  necessary”  disloca¬ 
tions  by  Ashby,  and  occur  in  order  to  allow  compatible  de¬ 
formation  of  a  system  with  geometrical  constraints  such  as 
hard  particles  which  do  not  deform  as  the  surrounding  duc¬ 
tile  matrix.  These  geometrically  necessary  dislocations  are 
required  if  the  deformation  takes  place  without  the  for¬ 
mation  of  voids  about  the  hard  particles.  Slip  dislocations 
are  a  function  of  the  material  properties  of  the  system,  and 
are  not  dependent  on  the  microstructural  constraints.  Ac¬ 
cording  to  Ashby,  the  density  of  geometrically  necessary 
dislocations,  p^.  is  given  by: 


12] 


where  A is  the  "geometrical  slip  distance”  analogous  to  the 
slip  distance  in  pure  crystals.  For  platelet  particles,  A*^  is 
equal  to  the  length  of  the  plate  and  y  is  the  shear  strain.  For 
a  1  pet  shear  strain  and  A”  =  4  pm.  equals  approxi¬ 
mately  3  X  10"  m  ^  Taking  these  dislocations  into  account 
results  in  a  further  addition  to  the  predicted  dislocation  den¬ 
sity  in  the  Al  matrix. 

The  purpose  of  this  investigation  was  to  determine  if 
dislocation  generation  occurs  at  the  Al/SiC  interface  on 
cooling  a  composite  from  annealing  temperatures,  and  to 
determine  if  the  observed  densities  of  dislocations  generated 
during  cooling  are  in  agreement  with  densities  predicted  by 
theoretical  calculations. 


II.  MATERIAL 

Three  types  of  material  were  examined.  The  first  was  a 
SiC  whisker  composite  purchased  from  ARCO  Silag.  It  is  a 
powder  metallurgy  product:  6061  Al  alloy  powder  is  com¬ 
pressed  with  SiC  whiskers  to  form  a  billet;  then  the  billet  is 
extruded  to  align  the  whiskers  and  form  a  12.5  mm  diameter 
rod.  Three  different  volume  fractions  of  SiC  were  consid¬ 
ered:  0  pet,  5  pet,  and  20  pet.  The  zero  vol  pet  material 
served  as  a  control. 

The  second  type  of  composite  was  purchased  from  DWA, 
and  contains  SiC  of  platelet  morphology.  The  platelets  are 
5  to  7  pm  long  and  have  an  aspect  ratio  of  two  to  three.  This 
composite  is  also  a  powder  metallurgy  product  supplied  in 
the  form  of  a  plate. 

Third  type  of  material  was  wrought  1100  grade  Al  in  the 
form  of  a  12.5  mm  diameter  rod,  and  it  was  in  the  as- 
received  condition.  This  material  also  served  as  a  control. 


ni.  SAMPLE  PREPARATION  AND 
EXAMINATION  PROCEDURE 

An  ion  milling  technique  w  as  required  for  the  production 
of  TEM  samples  due  to  the  SiC  in  the  Al  matrix. 

These  two  types  of  composite  and  the  0  vol  pet  control 
were  machined  into  rods  (12  mm  in  diameter,  4  cm  long), 
annealed  for  12  hours  at  a  solutionizing  temperature  of 
810  K,  and  furnace  cooled.  After  annealing,  slices  of 
0.76  mm  thickness  were  cut  by  electric  discharge  machin¬ 
ing  (EDM)  at  80  to  100  V.  Deformation  damage  from  EDM 
is  estimated  to  extend  0.20  mm  beneath  the  surface.'^  The 
slices  were  fixed  to  a  brass  block  with  double-sided  tape  and 
surrounded  by  brass  shims,  then  mechanically  thinned  on  a 
rotating  water  flooded  wheel  covered  with  400  then  600  grit 
paper  to  remove  the  EDM  damage  and  reduce  the  thickness 
to  approximately  0.127  mm.  Final  thinning  was  carried  out 
using  argon  ion  plasma  bombardment,  operating  at  6  kV, 
and  ion  current  of  50  micro  amperes  and  a  sample  inclina¬ 
tion  of  15  deg  to  the  ion  beam.  For  these  operating  parame¬ 
ters  the  projected  range,  or  average  distance  the  argon  ion 
travels  into  the  foil,  is  only  20  nm. ^  Dupuoy^'  conducted 
an  in  situ  ion  thinning  experiment  on  Fe  and  Al-Ag  speci¬ 
mens  using  a  3  MV  microscope.  Dislocation  arrangements 
and  microstructures  in  Al-Ag  and  Fe  are  not  altered  by  ion 
thinning  even  though  some  point  defects  are  introduced  into 
the  near  surface  region  of  the  sample  by  ion  bombardment. 
TTierefore,  it  can  be  concluded  that  ion-milling  does  not 
introduce  or  remove  dislocations  in  the  TEM  foils. 

The  1100  grade  Al  control  samples  were  prepared  from 
the  as-received  wrought  rod  in  the  same  manner  as  the 
composite  samples,  except  electro-polishing  was  employed 
instead  of  ion-thinning. 

The  thinned  samples  were  then  observed  in  the  HVEM 
operating  at  800  kV  with  a  beam  current  of  2.3  p.A.  A 
double  tilt,  side-entry,  furnace  type  heating  stage  was  used 
to  heat  the  specimen.  While  being  observed  in  the  micro¬ 
scope,  the  samples  were  heated  to  800  K  and  held  for  15 
minutes,  then  cooled  to  ambient  temperature.  Subgrains 
exhibiting  dislocations  in  contrast  were  chosen  for  obser¬ 
vation.  During  heating  and  cooling,  thermal  drift  of  the 
stage  and  thermal  expansion  and  contraction  of  the  sample 
caused  the  chosen  subgrain  to  move.  In  order  to  maintain  the 


same  subgrain  in  the  field  of  view  at  the  same  crystal¬ 
lographic  orientation,  it  was  necessary  to  slightly  translate 
and  tilt  the  specimen  almost  continuously.  Since  changing  to 
SAD  conditions  during  cycling  to  monitor  orientation  would 
have  resulted  in  loss  of  the  chosen  area  from  the  field  of 
view,  the  orientation  was  maintained  constant  by  monitoring 
the  contrast  of  microstructures  such  as  a  subgrain  boundary 
or  SiC/Al  interface  in  the  bright-field  mode.  One  thermal 
cycle  required  about  one  hour,  and  most  of  the  samples  were 
observed  throughout  several  thermal  cycles.  The  thickest 
regions  of  the  sample  penetrated  by  the  beam  were  chosen 
for  observation,  and  at  operating  voltage  of  800  KV  the 
beam  will  penetrate  0.8  thick  Al.^‘  The  dislocation  den¬ 
sity.  p,  was  determined  by  a  line  intercept  method  adapted 
from  Hale.^’  A  grid  of  lines  is  placed  over  the  TEM  micro¬ 
graph,  then  the  intersections  of  dislocation  lines  with  the 
grid  lines  are  counted.  The  dislocation  density,  p,  is  given 
by 


where  N  is  the  number  of  dislocation  intersections  with  the 
grid  lines,  L  is  the  length  of  the  grid  lines  divided  by  the 
magnification,  /  is  the  thickness  of  the  sample  (0.8  pm), 
and  the  length  of  the  lines  on  the  grids  was  0.58  m.  Each 
reported  dislocation  density  is  an  average  value  obtained 
from  3  to  10  micrographs. 


IV.  RESULTS 

The  discussion  of  the  experimental  results  will  be  divided 
into  three  parts;  0  vol  pet,  5  vol  pet,  and  20  vol  pet  of  SiC. 
In  all  cases,  typical  results  will  be  described.  A  total  of  800 
micrographs  were  taken. 

A.  Controls.  0  Vol  Pet  SiC.  6061  Al,  and  1100  Al 

The  0  vol  pet  6061  Al  control  sample  had  a  large  sub¬ 
grain  size  (approximately  5  pm)  and  a  low  dislocation  den¬ 
sity  (8  X  10'^  m’^)  (Figure  1(a)).  A  few  of  the  larger 
second  phase  precipitate  particles  (Mg;Si)  were  surrounded 
by  dislocation  tangles  bowing  out  from  the  precipitate  inter¬ 
face  (Figure  1(b)).  On  heating,  the  dislocations  began  to 
move  in  the  sample,  migrating  away  from  the  particles 
Other  dislocations  were  also  generated  in  other  areas  of  the 
subgrain  and  moved  through  the  matrix,  occasionally  being 
pinned  by  precipitate  panicles.  Eventually,  at  elevated  tem¬ 
peratures  (670  K),  all  of  the  dislocations  disappeared.  Slip 
traces  left  behind  when  the  dislocations  moved  also  disap¬ 
peared  at  temperatures  close  to  700  K,  possibly  due  to  sur¬ 
face  diffusion.  The  sample  was  held  at  800  K  for  about 
15  minutes. 

The  image  tended  to  be  out  of  focus  at  high  temperatures 
due  to  thermal  drift  of  the  stage  (Figures  1(b)  and  1(c)).  On 
cooling,  dislocations  reappeared  at  the  large  particle  inter¬ 
face  at  about  500  K,  sometimes  moving  faster  than  they 
could  be  seen.  The  dislocations  formed  tangles  in  the  vicin¬ 
ity  of  the  precipitate;  however,  most  of  the  matrix  did  not 
accumulate  any  dislocations  (Figure  1(d)).  There  were  a 
few  dislocations  generated  at  subgrain  or  grain  boundaries. 

The  1 100  Al  sample  had  a  dislocation  density  initially  of 
4  X  10'^  m and  a  large  subgrain  size,  5  |im  (Figures  2(a) 


and  2(b)).  Thicker  regions  containing  few  precipitates  were 
chosen  for  observation.  On  heating,  most  of  the  dislocations 
had  disappeared  by  673  K.  At  this  temperature,  the  heating 
stage  mechanism  failed;  therefore,  the  sample  never  reached 
800  K,  but  heating  was  sufficient  to  remove  the  dislocations 
from  the  area  under  observation.  The  dislocations  in  this 
area  did  not  return  on  cooling,  except  a  few  which  were 
connected  to  precipitate  particles  (Figures  3(a)  and  3(b)). 

B.  5  Vol  Pet — Transverse  and  Longitudinal  Whisker  SiC 

The  5  vol  pet  whisker  sample,  sectioned  transverse  to  the 
whiskers  so  that  the  hexagonal  whiskers  are  viewed  end-on 
in  the  microscope,  had  a  small  subgrain  size,  2.0  /urn,  and 
a  high  dislocation  density.  4  x  10”  m“^  in  the  subgrains 
(Figure  4(a)).  In  this  case,  on  heating  the  dislocations  did 
not  disappear  but  straightened  from  an  initially  bowed  con¬ 
figuration  and  became  more  regularly  spaced;  a  polygonized 
configuration.  These  dislocations  did  not  disappear  at  high 
temperatures  (Figures  4(b)  and  4(c)).  On  cooling,  the  dis¬ 
locations  again  bowed  away  from  the  Al/SiC  interface  but 
the  density  of  dislocations  did  not  increase;  if  anything,  the 
number  slightly  decreased  (Figures  4(d)  and  4(e)).  Several 
of  the  subgrains  appeared  to  change  shape  during  the  ther¬ 
mal  cycle;  however,  their  boundaries  did  not  move  past  the 
surrounding  SiC  whiskers. 

The  longitudinal  5  vol  pet  sample,  sectioned  parallel  to 
the  SiC  whisker  axis  so  that  the  whiskers  could  be  viewed 
lengthwise,  did  show  the  characteristic  disappearance  of  dis¬ 
locations  at  high  temperatures,  and  then  the  return  of  dislo¬ 
cations  on  cooling,  especially  at  the  ends  of  the  whiskers 
(Figures  5(a)  through  (f)).  The  parallel  lines  in  the  whiskers 
have  been  previously  identified  as  microtwins.” 

C.  20  Vol  Pet  Whisker  and  Platelet 

The  microstructure  of  the  20  vol  pet  whisker  sample  be¬ 
fore  heating  is  characterized  by  a  small  subgrain  size  of  the 
order  of  2  to  3  /zm  and  also  by  a  dislocation  density  of  about 
10”  m"^  (Figure  6(a)).  The  sample  was  heated  while  focus¬ 
ing  on  a  single  subgrain  surrounded  by  several  SiC  whisk¬ 
ers.  The  dislocations  began  to  move  and  rearrange,  some 
moving  very  quickly,  and  eventually  disappearing.  Upon 
reaching  470  K  most  of  the  dislocations  had  disappeared. 
The  sample  was  heated  further  to  800  K  (Figure  6(b)),  held 
fora  few  minutes,  and  then  cooled.  On  cooling,  dislocations 
reappeared  and  emanated  from  the  whisker  interface  form¬ 
ing  a  tangle  of  dislocations  in  the  small  subgrain.  Although 
some  dislocations  began  appearing  on  cooling  at  about 
473  K,  a  great  number  of  them  formed  at  temperatures  less 
than  373  K.  The  dislocation  density  after  cooling  was  com¬ 
parable  to  the  dislocation  density  before  the  thermal  cycle 
(Figures  6(c),  6(d),  and  6(e)).  When  black  spots  began  ob¬ 
scuring  the  picture,  another  area  free  of  the  black  spots  was 
moved  into  the  field  of  view  . 

The  20  vol  pet  platelet  sample  (Figure  7(a))  exhibited  be¬ 
havior  similar  to  the  whisker  sample,  in  that  most  of  the 
dislocations  disappeared  upon  reaching  650  K  (Figure  7(b)). 
and  then  on  cooling,  dislocations  reappeared  (Figures  7(c) 
and  7(d)).  However,  certain  aspects  of  the  microstructural 
alterations  on  cooling  were  peculiar  to  the  platelet  sample. 
In  the  first  cycle,  the  dislocations  disappeared  and  then 
reappeared  on  cooling.  An  unusual  subgrain  was  observed 
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Fig  I  —  (a)  Micrograph  of  a  control  sample  illustrating  the  low  dislocation  densit>'  and  the  large  grain  size  of  annealed  powder  compacted  6061  Al  with 
no  SiC.  The  few  dislocations  present  in  this  control  sample  are  associated  with  the  second  phase  precipitates  of  the  6061-AI  alloy  Compare  the  dislocation 
densities  and  subgrain  sizes  with  those  of  composites  shown  in  Figs.  4(a),  6(a),  and  7(a).  (b)  A  large  second  phase  precipitate  (top  center)  surrounded  by 
dislocations  at  the  beginning  of  the  thermal  cycle  (340  K).  (c)  As  temperature  increased  to  annealing  temperatures  in  the  microscojie  (780  K),  the  dislocations 
gradually  disappeared  until  only  a  few  remained  (same  area  as  Fig  1(b)).  (</)  During  in  situ  cooling  of  the  control  sample  dislocations  reformed  in  the  area 
of  the  precipitates,  shown  here.  However,  most  of  the  matrix  remained  free  of  dislocations  (same  area  as  Figs  1(a)  and  1(b)). 


which  was  filled  with  slip  traces.  This  area  was  chosen  as 
the  area  of  focus  in  the  next  cycle.  The  parallel  lines  disap¬ 
peared  on  heating;  then  on  subsequent  cooling,  packets  of 
slip  traces  appeared  emanating  from  the  SiC  platelet.  The 
slip  traces  formed  in  three  directions  at  angles  of  about 
82  deg  and  45  deg  to  each  other  and  were  associated  with 
dislocations  which  had  appeared.  The  dislocations  causing 
the  slip  traces  had  moveti  too  quickly  to  be  seen.  After  the 
third  thermal  cycle,  many  subgrains  were  filled  with  slip 
traces  originating  at  the  SiC  platelet  interface, 

D,  Dislocation  Densities 

The  dislocation  densities  of  all  samples  before  and  after 
in  situ  annealing  are  shown  in  Table  1,  The  densities  are 
higher  in  tfie  composite  samples  than  in  the  control  samples 


both  before  and  after  in  situ  thermal  cycling,  i.e.,  an¬ 
nealing,  However,  the  dislocation  densities  reported  are 
lower  limit  densities.  There  was  difficulty  in  taking  selected 
area  diffraction  patterns  of  a  given  subgrain  and  then  tilting 
to  a  specific  reflection,  i.e.,  L420J  and  assuring  that  the 
same  orientation  was  maintained  throughout.  Therefore,  the 
reported  densities  could  be  Vi  to  'A  of  the  actual  density.*' 
The  distribution  of  the  dislocations  within  the  samples 
was  not  uniform;  there  was  a  higher  density,  as  to  be  ex¬ 
pected,  near  the  SiC.  Also  greater  dislocation  generation 
was  observed  at  larger  SiC  platelets  and  whiskers  than  at 
the  smaller  ones,  and  more  dislocations  were  generated 
at  the  ends  of  the  SiC  whiskers  where  plastic  strain  dur¬ 
ing  cooling  is  greatest,  than  at  the  middle  of  the  whisker 
length.  The  dislocation  g  •deration  in  the  platelet  sample 
was  much  greater  than  the  (hslocation  generation  in  the 
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Fig  2  —  (a)  Wrought  1100  Al  control  sample,  as-received  condition.  Be¬ 
fore  in  situ  thermal  cycling,  this  control  sample  with  no  SiC  had  a  high 
dislocation  density  Like  the  6061  Al  control,  this  sample  has  a  large  grain 
iiae.  and  the  dislocations  present  are  often  associated  with  precipitates.  The 
white  structure  is  a  corrosion  tunnel  from  electropolishing  (h)  Low  mag¬ 
nification  view  of  1 100  control  sample  including  area  of  Fig.  2(a) 
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Fig.  3  —  (a)  The  same  areas  in  Fig  2(a)  after  in  situ  thermal  cycling  In 
the  absence  of  SiC,  few  dislocations  reappear  in  a  thermally  cycled  sample, 
the  few  dislocations  present  arc  associated  with  second  phase  precipi¬ 
tate  panicles. (fc)  Low  magnification  view  of  area  in  Fig  3(a)  after  ther¬ 
mal  cycling 
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Fig.  4  —  (o)  A  foil  cut  transverse  to  the  whisker  axis  (5  vol  pci  SiC).  At 
the  beginning  of  the  thermal  cycle.  370  K.  the  dislocations  tended  to  bow 
away  from  the  two  hexagonal  SiC  whiskers  (fc)  The  dislocations  straighten 
on  further  heating  to  570  K  (same  area  as  Fig.  4(a)).  (r)  Further  re¬ 
arrangement  of  the  dislocations  occurred  on  heating  to  the  armealing  tem- 
peramre  of  800  K  in  the  nucroscope.  The  dislocations  straightened  and 
moved  to  more  equally  spaced  positions  (same  area  as  Figs.  4(a)  and  4(b)l. 
(<f)  During  cooling  to  570  K  the  dislocations  again  began  to  bow  away 
from  the  SiC  whiskers  (same  area  as  Figs  4(a-c))  (e)  On  cooling  to 
360  K.  the  curvature  of  the  dislocations  became  pronounced,  indicating  the 
operation  of  a  thermal  stress;  however,  the  generation  of  new  dislocations 
did  not  occur  (same  area  as  Figs  4(a-dl). 


whisker  samples,  for  SiC  platelets  are  about  5  times  larger 
than  SiC  whiskers. 

As  a  consequence,  of  all  the  dislocation  generation,  there 
is  a  possibility  of  void  formation  at  the  SiC-Al  interface. 
Voids  were  not  observed  in  any  of  the  samples. 


DISCUSSION 


The  presence  of  SiC  particles  of  either  whisker  or  platelet 
morphology  in  an  A1  metal  matrix  composite  resulted  in  the 
generation  of  dislocations  at  the  Al/SiC  interface  when  the 
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Fig  5  —  (a)  This  sample  was  prepared  from  the  same  5  vol  pci  whisker 
material  shown  in  Figs  4(a-e).  only  sectioned  parallel  to  the  SiC  whisker 
axis  Dislocations  are  associated  with  the  SiC  and  second  phase  precipitates 
at  the  beginning  of  the  thermal  cycle  at  340  K  Also  visible  are  parallel  slip 
traces  from  the  previous  thermal  cycle  (b)  The  same  region  at  annealing 
temperature.  810  K.  after  most  dislocations  have  disappeared  (c)  The 
same  area  after  one  dIsliKation  has  reappeared  on  cooling  to  580  K  .  (d )  On 
further  cooling  the  sample  to  480  K  more  dislocations  have  formed,  some 
leaving  slip  traces  as  they  moved  away  from  the  Al/SiC  interface,  (e)  On 
cooling  to  430  K.  dislocation  activity  greatly  increased  as  evidenced  by  the 
higher  density  of  dislocations  and  slip  traces,  (f)  Another  location  in  the 
sample  at  340  K  which  also  shows  dislocation  generation  clearly  associated 
with  the  SiC.  This  area  was  not  irradiated  by  the  electron  beam  during  the 
thermal  cycle. 
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composite  was  cooled  from  the  annealing  temperature.  In 
general,  the  high  density  of  dislocations  originally  present  in 
the  composite  samples  disappeared  at  500  to  650  K,  then 
reappeared  on  cooling  at  densities  close  to  the  high  densities 
originally  observed  in  the  annealed  specimens.  In  contrast, 


id) 

Fig  6 — (a)  The  iniiial  high  dislocation  density  and  small  subgrain  sire  of 
the  20  vol  pci  whisker  composites  is  shown  here  (b)  The  dislocations 
disappeared  from  the  sample  dunng  in  situ  heating  to  800  K  Thermal  drift 
of  the  stage  prevented  sharp  focusing  of  the  image,  (r)  On  cooling  to 
37S  K.  dislocations  had  reappeared  in  the  sample,  emanating  from  the 
AI/SiC  interface  and  moving  into  the  matrix  until  a  dislocation  tangle 
formed,  (d)  The  clear  area  at  375  K  (Fig  6(c))  as  temperature  decreased 
to  340  K  Black  spots  also  began  to  form  on  the  sample  at  this  tempera¬ 
ture.  (e)  At  lower  center,  the  appearance  of  a  hook-shaped  dislocation  is 
observed.  The  black  spots  have  become  more  prominent.  330  K 


dislocations  were  not  generated  to  the  same  extent  in  either 
of  the  control  samples  during  cooling. 

The  small  subgrain  size  and  high  dislocation  densities 
previously  observed  by  Arsenault  and  Fisher’  can  be  associ¬ 
ated  with  the  presence  of  the  SiC  in  the  composite,  specifi¬ 
cally  the  CTE  differential .  The  difference  in  CTE  resulted  in 
stresses  large  enough  to  cause  plastic  deformation,  i.e..  the 
generation  of  dislocations.  These  dislocations  can  be  de¬ 
fined  as  slip  dislocations.  Dislocation  generation  is  also 
required  to  accommodate  the  heterogeneous  plastic  flow 
in  the  vicinity  of  the  deforming  matrix  since  voids  are 
not  observed.  We  can  conclude,  therefore,  that  both  of 
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Fig  7  —  (u  I  The  diffraction  contrast  of  the  high  dislocation  density  in  the  central  subgrain  darkens  the  subgrain  in  this  20  vol  pci  SiC  platelet  composite. 
The  subgrain  is  surrounded  by  at  least  three  irregularly  shaped  SiC  platelets  which  are  approximately  the  same  size  as  the  subgrain,  lb)  During  heating  to 
6^0  K  most  of  the  dislocations  disappeared  from  the  sample.  Ihe  pebble-like  second  phase  precipitates  can  be  clearly  distinguished  here.  The  same  area  is 
shown  in  Id)  filled  with  slip  traces  (r)  On  cooling  to  340  K.  dislocations  reappeared,  resulting  in  the  formation  of  dense  slip  traces  showing  the  paths  of 
the  disliK'ations  as  they  moved  across  or  out  of  the  sample.  (</)  This  is  the  same  area  as  shown  in  (b)  after  cooling  to  340  K  The  area  was  free  of  almost 
all  disliKalions  at  high  temperature  The  subgrain  is  now  filled  with  slip  traces  and  disliKations  The  slip  traces  are  clearly  associated  with  the  SiC  platelets 
indicating  that  the  ongin  of  the  dislocations  seen  in  the  sample  at  low  temperature  is  the  Al/SiC  interface 


Table  I.  Dislocation  Density  before  and  after  Thermal  Cycling  of  Samples;  Dislocation  Density  (m'‘) 


Sample 

Before  Thermal  Cycle 

After  Thermal  Cycle 

( 1 )  20  vol  pet  whisker* 

2.0  X  10" 

1.0  X  10" 

(2)  20  vol  pci  platelet* 

2.0  X  10" 

4.0  X  10'- 

(3|  5  vol  pet  whisker* 

4.0  X  10" 

2.0  X  10" 

(4)  0  vol  pet  whisker’ 

8  0  X  lO'- 

6.4  X  10'- 

(5)  1100  Al  sample* 

4.0  X  10'- 

5.3  X  10" 

'Bulk  annealed  12  hours  at  800  K 
•As-received  1100  grade  wrought  Al  alloy 

these  mechanisms  are  operating,  which  results  in  a  high  found.  A  stable  polygonized  substructure  appeared  to  pre¬ 
dislocation  density  due  to  the  SiC  in  the  matrix.  In  the  vent  the  disappearance  of  the  dislocations  at  high  tern- 

vicinity  of  precipitates,  transformation  strains  can  also  pro-  peratures.  Backstresses  from  these  dislocations  on  cooling 
duce  dislocations.  could  have  prevented  further  generation  of  dislocations  in 

The  exception  to  the  generally  observed  disappearance  the  small  subgrains  on  cooling, 
of  dislocations  on  heating  was  the  5  vol  pet  transverse  The  intensity  of  dislocation  generation  at  the  SiC- A1  inter¬ 

sample  for  which  no  satisfactory  explanation  has  been  face  is  related  to  size  and  shape  of  SiC  particle .  The  intensity 


of  generation  is  lowest  for  small,  nearly  spherical  particle. 
As  the  particle  size  increases,  i.e.,  from  1  /xm  to  5  fim,  the 
intensi^  of  generation  increases  significantly.  Also  the  in¬ 
tensity  is  much  greater  at  the  comer  of  a  particle  than  along 
the  sides  as  is  evidenced  in  the  longitudinal  whisker  and 
platelet  composite  samples.  The  relationship  between  the 
size  of  the  particle  and  the  plastic  zone  has  been  qual¬ 
itatively  described  by  Lee  et  al.'^  and  they  predicted  that  as 
the  particle  size  increases  the  plastic  zone  size  increases. 
Also,  Lee^‘  has  shown  that  the  plastic  strain  about  the  comer 
of  a  particle  is  greater  than  along  the  side  (which  should  be 
intuitively  obvious). 

Experimental  conditions  which  may  influence  the  results 
must  be  recognized  and  considered  in  an  HVEM  experi¬ 
ment,  because  the  sample  is  exposed  to  high  energy  elec¬ 
trons,  and  thin  foil  samples  are  used  to  approximate  bulk 
behavior.  Surface  effects  are  among  the  most  important 
effects  to  be  considered,  since  it  is  easy  for  dislocations 
generated  at  the  Al/SiC  interface  to  move  out  of  the  sample 
through  the  surfaces  of  a  TEM  foil.  Thick  sections  of  the  foil 
were  examined  in  the  HVEM  in  order  to  reduce  surface 
effects;  nevertheless,  dislocation  relaxation  out  the  surface 
occurred,  and  resulted  in  a  reduction  of  the  observed  dis¬ 
location  density.  This  effect  is  most  apparent  in  the 
20  vol  pet  platelet  sample  where  the  slip  traces,  indicating 
that  dislocations  have  moved  out  of  the  sample,  nearly  cover 
the  entire  surface  of  the  sample.  In  bulk  specimens,  dis¬ 
locations  would  accumulate  in  the  subgrains  until  the  en¬ 
suing  backstresses  due  to  the  pile-up  exceeded  the  local 
yield  stress  surrounding  the  particle.  Also,  the  geometry  of 
a  thin  foil  specimen  allows  elastic  relaxation  of  stresses  on 
cooling  by  buckling,  also  giving  an  artificially  low  value  of 
the  dislocation  density. 

The  effects  of  irradiation  of  the  samples  by  the  high 
energy  electrons  of  the  beam  must  also  be  considered.  Elec¬ 
tron  irradiation  of  the  sample  can  result  in  the  formation  of 
vacancy  clusters  and  small  dislocation  loops  which  appear 
as  black  spots  and  then  grow  to  form  dislocation  tangles." 
The  observation  of  black  spots  on  some  of  the  samples 
(Figures  1(d)  and  6(d))  indicated  that  electron  irradiation 
damage  most  probably  occurred. 

The  control  samples  were  invaluable  in  determining  that 
the  dislocations  generated  on  cooling  were  not  artifacts  due 
to  the  effects  of  electron  irradiation.  In  the  6061  0  vol  pet 
SiC  control  sample  a  few  dislocations  were  generated  at  a 
few  large  precipitates. 

The  1  lOO  control  samples  were  also  exposed  to  the  high 
energy  electron  beam  and  substantial  formation  of  dis¬ 
locations  did  not  occur  (Figures  3(a)  and  3(b)).  Also  the 
intensity  of  dislocation  generation  can  be  correlated  with  the 
size,  volume  fraction,  and  shape  of  the  SiC  or  second  phase 
precipitate  particles  present,  indicating  that  the  particles  and 
not  the  electrons  of  the  beam  were  the  cause  of  the  dis¬ 
locations  being  generated.  A  more  likely  explanation  of  the 
appearance  of  the  black  spots  is  beam  contamination.  The 
decontaminator  was  not  always  operating,  and  beam  con¬ 
tamination  usually  condenses  on  samples  at  temperatures 
less  than  473  K  which  coincides  with  our  observations.^’-^* 
Some  of  the  black  spots  could  also  be  due  to  second  phase 
precipitation  in  the  samples  on  cooling,  since  the  spots  were 
sometimes  preferentially  associated  with  interfaces  and 
grain  boundaries  (Figures  1(d)  and  7(d)). 


Beam  heating  is  another  factor  to  be  considered.  For  our 
operating  conditions,  beam  heating  of  the  sample  is  approxi¬ 
mately  10  to  IS  K,  and  this  could  have  had  an  effect  on 
dislocation  generation  in  the  sample  due  to  the  thermal  gra¬ 
dient  which  is  induced.’”^  But  examination  of  the  same 
composite  samples  without  thermal  cycling  did  not  result  in 
dislocation  generation  at  the  Al/Si(r  interface,  and  these 
samples  were  exposed  to  the  same  beam  conditions. 

Due  to  the  difficulty  associated  with  tilting  a  very  fine 
subgrained  material  to  the  various  diffracting  conditions 
required  to  image  all  of  the  dislocations  in  the  subgrain,  the 
reported  densities  could  be  Vy  to  '/:  below  the  actual  number. 
Although  more  rigorous  tilting  would  give  more  precise 
values,  a  good  idea  of  the  relative  densities  in  the  samples 
can  be  obtained  by  imaging  dislocations  in  many  subgrains 
for  each  material  and  assuming  that  the  value  will  be  sys¬ 
tematically  low  for  all  the  samples. 

The  net  result  of  considering  all  of  the  experimental  fac¬ 
tors  which  may  influence  the  experimentally  determined 
dislocation  densities  after  a  thermal  cycle,  is:  (1)  the  dis¬ 
location  generation  observed  during  cooling  can  be  readily 
attributed  to  differential  thermal  contraction  of  the  A1  and 
SiC,  (2)  the  observed  densities  are  lower  than  the  densities 
which  would  be  observed  if  bulk  samples  could  be  exam¬ 
ined  and  if  diffracting  conditions  were  controlled  to  image 
all  dislocations. 

It  should  also  be  pointed  out  that  slip  line  generation 
about  a  SiC  cylinder  in  an  A1  disk  due  to  thermal  cycling  has 
been  demonstrated  by  Flom  and  Arsenault.’® 

VI.  CONCLUSIONS 

From  a  consideration  of  the  experimental  results,  the 
following  conclusions  can  be  drawn: 

1.  The  high  dislocation  density  (10"  m"’)  previously  ob¬ 
served  in  bulk  annealed  composites  is  due  to  differential 
thermal  contraction  of  Al  and  SiC  on  cooling  from  the 
elevated  temperatures  of  annealing. 

2.  The  density  of  dislocations  observed  in  this  experiment 
as  a  result  of  thermal  cycling  is  lower  than  the  actual 
density  generated  during  thermal  cycling  because  dis¬ 
locations  are  lost  through  the  surfaces  of  the  thin  foil 
samples  during  cooling. 

3.  The  densities  of  dislocations  observed  (10"  m"’)  would 
be  equal  to  the  high  densities  previously  observed  in 
bulk  annealed  composites  if  it  were  not  for  dislocation 
loss  through  the  surfaces,  and  the  observed  densities 
would  be  closer  to  densities  predicted  by  calculations 
♦  X  10'*  m"’,  if  it  were  not  for  dislocation  loss  through 
the  surfaces. 

4.  Thermal  cycling  causes  the  disappearance  of  dislocations 
at  high  temperatures  and  the  generation  of  dislocations  at 
Al/SiC  interfaces  and  precipitates  on  cooling.  Subgrain 
growth  is  hindered  by  the  presence  of  SiC  particles. 
Transformation  strains  also  cause  dislocation  generation 
at  precipitates,  and  polygonized  configurations  prevent 
the  disappearance  of  dislocations  at  high  temperatures. 
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Abstract 


When  a  metal  matrix  composite  Is  cooled  down  to  room 
temperature  from  the  fabrication  or  annealing  temperature, 
residual  stresses  are  Induced  in  the  composite  due  to  the 
mismatch  of  the  thermal  expansion  coefficients  between  the 
matrix  and  fiber. 

An  Investigation  was  undertaken  of  the  magnitude  of  the 
thermal  residual  stresses  by  determining  the  difference  of 
the  yield  stresses  (Ao  )  between  tension  and  compression 
resulting  from  the  thermal  residual  stresses.  A  theoretical 
model  based  on  the  Eshelby’s  method  was  then  constructed  for 
the  prediction  of  the  thermal  residual  stresses  and  Aa  . 

The  agreement  obtained  was  very  good  between  the  experimental 
results  and  the  theoretical  predictions. 


"n  ^  iji  gi  ijii  ui'j  tP'>t"*y 


»L»-r '-  w  •■»■ 


2 


1.  Introduction 


! 


Metal  matrix  composites  (MMCs),  Including  eutectic 
composites,  are  becoming  important  In  their  application  to 
structural  components  which  are  to  be  used  at  intermediate 
and  high  temperatures.  When  MMCs  are  fabricated  at  high 
temperature  or  annealed  at  a  certain  high  temperature  and 
cooled  down  to  room  temperature.  It  results  In  undesirable 
properties,  such  as  low  tensile  yield  and  ultimate  strengths. 
These  results  are  mainly  due  to  residual  stresses  that  are 
car.sed  by  the  mismatch  of  the  thermal  expansion  coefficients 
between  the  matrix  and  fiber.  This  subject  has  been  studied 
by  a  number  of  researchers  (for  examples  see  references  1- 
13).  The  residual  stresses  so  induced  have  been  observed  In 
tungsten  fiber/copper  composites and  In  SIC  whlsker/606l 
A1  composltes^®*^^.  Most  of  the  models  purposed  to  estimate 
the  residual  stress  were  based  on  one-dlmenslonal 
analysis^ continuous  fiber  system^^»^\  or  spherical 
particle  system^^®^. 

The  model  Is  based  on  Eshelby's  equivalent  Inclusion 
method^'^\  which  has  been  used  to  solve  the  problem  of 
thermal  residual  stress ^ \  The  advantages  of  the 
Eshelby's  model  are  that  It  can  solve  a  three-dimensional 
composite  system  such  as  short  whisker  composite  and  also  can 
take  Into  account  the  effect  of  the  volume  fraction  of  fiber 
(Vj.)  easily. 
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Eshelby's  method  has  been  also  used  for  predicting  the 
yield  stress  and  work-hardening  rate  of  metal  matrix 
composites^^^”^^^ .  The  effect  of  the  thermally  Induced 
residual  stress  on  the  yield  stresses  has  been  discussed  by 
Wakashlma  et  al.^^®^  who  predicted  that  the  yield  stress  in 
compression  (a^)  exceeds  that  In  tension  (o^)  for  continuous 

y  y 

B  flber/Al  composites  which  were  cooled  down,  although  no 
comparison  with  the  experimental  data  was  made. 

In  this  paper  we  focus  on  the  residual  stresses  Induced 
In  a  short  whisker  MMC  due  to  the  temperature  drop  and  its 
effect  on  the  yield  stresses.  The  target  MMC  is  aligned 
fiber  MMC  with  a  special  emphasis  on  SIC  whlsker/606l  A1 
composite. 


L’ 


■’*■1*1.**  ■  1*  *1^.  M 
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2.  Experimental  Procedure 

A.  Materials 

The  target  short  fiber  MMC  was  SIC  whlsker/606l  A1 
composite  and  was  purchased  from  ARCO/SILAG.  As  a  control 
MMC  a  SIC  spherical  1100  A1  composite  was  also  purchased  from 
ARCO/SILAG.  The  composite  was  In  the  form  of  an  extruded  rod 
15.5  mm  In  diameter.  Three  different  volume  fractions  of 
fiber  (or  whisker)  composites  were  used:  =  0,  0.05  and 

0.2,  and  a  composite  with  spherical  SIC  particles  (0.5 
In  diameter  having  a  volume  fraction,  Vg,  equal  to  0.2  was 
also  used.  All  of  the  material  was  supplied  as-fabricated 
(no  heat  treatment)  condition,  they  then  were  machined  into 
samples  (Fig.  la  and  lb),  annealed  for  12  hours  at  810  K,  and 
then  furnace  cooled. 

B.  Testing  Methods 

The  tension  and  compression  tests  were  performed  in  an 
Instron  testing  machine  using  a  liquid  metal  container  as  the 
lower  gripping  device.  This  method  of  gripping  was  employed 
to  ensure  very  good  alignment  of  the  sample.  This  test 
procedure  Is  described  In  greater  detail  elsewhere^^®^ . 
Samples  of  two  different  gauge  lengths  were  used  to  determine 
If  there  was  a  gauge  length  effect  (there  was  none).  The 
effective  gauge  length  of  the  samples  was  determined  by 
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glueing  strain  gauges  on  the  center  portion  of  the  sample  and 

comparing  the  results  obtained  from  a  clip  on  extensoraeter 

which  was  mounted  Into  the  "V”  grooves.  If  the  extensoraeter 

is  mounted  directly  to  the  uniform  gauge  section,  there  is  a 

high  probability  that  the  sample  will  fracture  where  the 

"knife"  edge  of  the  extensoraeter  makes  contact  with  the 

sample.  Several  tests  in  the  low  stress  range  were  conducted 

using  both  the  extensoraeter  and  the  strain  gauge;  from  these 

tests  the  effective  gauge  length  was  determined.  Subsequent 

tests  were  conducted  using  only  the  extensometer .  The 

samples,  which  were  tested  in  the  range  from  8  x  10~^  sec”^ 

-3  -1 

to  2  X  10  sec  ,  showed  no  effect  on  strain  rate. 

The  conventional  X-ray  technlque^^^ ^  was  used  first,  but 
it  was  determined  that  there  were  problems  when  this  tech¬ 
nique  was  applied  to  discontinuous  composites.  The  residual 
stress  was  determined  by  measuring  the  change  In  lattice 
parameters.  There  Is  an  Implicit  assumption  made,  l.e.,  the 
matrix  remains  cubic. 


3.  Experimental  Results 


The  stress-strain  curves  of  the  Vj.  =  0,  0.05  and  0.2 
composites  are  plotted  In  Fig.  2a,  2b,  and  2c,  respectively. 
The  solid  and  dashed  curves  denote  the  tensile  and  compres¬ 
sive  test  results,  respectively.  The  yield  stress  was 
measured  as  0,2%  off  set  strain  and  Is  Indicated  In  the 
figures  by  arrows. 

The  shapes  of  the  stress-strain  curves  of  the  Vg  *  0.2 
composites  are  very  similar  to  those  of  whisker  composites, 
except  that  the  tension  stress-strain  curve  Is  below  that  of 
the  compression  stress-strain  curve.  However,  as  In  the  case 
of  Vj.  =  0,  there  was  no  difference  In  the  stress-strain 
curves  for  tension  and  compression  of  Vg  *  0. 

The  data  obtained  from  the  X-ray  analysis  is  given  In 
Table  I,  and  it  was  observed  that  there  Is  a  small  residual 
tensile  stress  In  the  matrix. 
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Theoretical  Procedure 

The  theoretical  model  used  here  Is  based  on  Eshelby's 
equivalent  Inclusion  model.  Mori  and  his  co-workers^ 
extended  Eshelby’s  method  to  predict  the  yield  stress  (a  ) 

y 

and  work-hardening  rate  of  aligned  short  whisker  composites. 
Wakashlma  et  al.^^®^  extended  the  above  approach  to  predict 
a  by  considering  the  mismatch  of  the  thermal  expansion 

y 

coefficients  of  the  matrix  and  fiber.  Following  the  above 
models,  Takao  and  Taya^^^^  have  recently  computed  the  stress 
field  In  and  around  a  short  fiber  In  a  short  fiber  composite 
where  the  fiber  Is  anisotropic  both  In  stiffness  and  thermal 
expansion. 

In  this  paper  we  focus  on  the  average  thermal  residual 
stress  Induced  In  the  matrix  by  the  cool-down  process  and 
also  the  yield  stresses  In  tension  (a^)  and  compression  (o*^) 

y  y 

when  the  composite  Is  tested  at  the  room  temperature.  The 
former  case  Is  essentially  based  on  the  model  by  Takao  and 
Taya^^^\  and  as  the  latter  case,  we  extend  the  model  by 
Wakashlma  et  al.^^®^  to  account  for  the  bl-llnear  stress- 
strain  curve  of  the  matrix. 


A.  Formulation 


Consider  an  Infinite  body  (D)  which  contains  ellipsoidal 
whiskers  (n )  aligned  along  the  x^-axls  (Pig.  3).  This  com¬ 
posite  body  D  Is  subjected  to  the  applied  stress  field  o?.. 
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The  stiffness  tensors  of  the  matrix  (D-n)  and  fiber  (ft)  are 
denoted  by  and  respectively.  Pollovrlng  Eshelby, 

the  transformation  straln^^^^  or  elgenstraln^^^^  Is  given  in 

K  • 

the  fiber  domain  U  as  where  is  the  strain  due  to  the 

mismatch  of  the  thermal  expansion  coefficients  and  the  uni¬ 
form  plastic  strain  e^^  is  prescribed  in  the  matrlx^^^^.  As 
far  as  the  stress  field  is  concerned,  the  model  of  Fig.  U  Is 
equivalent  to  that  of  Fig.  3.  Thus,  the  present  problem  Is 
reduced  to  "inhomogeneous  inclusions  problem"  (Pig. 

The  model  of  Fig.  ^  will  be  used  not  only  to  predict  the 
yield  stresses  and  work-hardening  rates,  but  also  to  compute 

the  thermal  residual  stresses.  For  the  latter  case,  we  will 

..  o  p 

set  0, .  «  ef.  “  0. 

IJ  IJ 

Following  the  Eshelby’ s  equivalent  inclusion  method 
modified  for  a  finite  volume  fraction  of  flbers^^^» the 
total  stress  field  in  the  fibers  is  given  by 


°1J  °1J  *  ^ljk£^®k£  ®k£  ®k£  "  ^®k£  ”  ®k£^^ 


_  f  O  ^  ^  D 

ie,  +  e,  e,  •  (a,  •  er 

ljk£^  k£  k£  k£  ^  kl  k£ 


where 


o  „  o 
®1J  *  ^ljk£®k£ 


in  D 


'^°lj^m  “  ^ljk£®k£ 


in  D-q 
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In  Eq.  (1)»  Is  the  fictitious  elgenstrain^ ^3)  which  was 

Introduced  to  connect  the  present  problem  to  "inclusion 
problem"  and  e^^j  Is  the  average  strain  disturbance  In  the 
matrix  and  related  to  the  average  stress  disturbance  In  the 
martlx  <o..>  by  Eq.  (3).  <o-.>  Is  defined  by 


where  Is  the  total  volume  of  the  matrix  and  and  ej^j 

In  Eqs.  (1)  and  (^)  Is  the  stress  and  strain  disturbance, 
respectively,  by  a  single  fiber  a  when  It  Is  embedded  In  an 
Infinite  matrix.  Since  the  stress  disturbance  when  It  Is 
Integrated  in  the  entire  domain  D  vanishes. 


(8) 
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According  to  Eshelby,  Is  related  to  the  total 


elgenstraln  ej^^  as 


•  « 


where  Is  the  Eshelby 's  tensor  and  a  function  of 


^IJki  geometry  of  the  ellipsoidal  flber^^^^.  The 


fiber  is  assumed  as  a  prolate  spheroid,  hence  the  fiber 


aspect  ratio  t/d  is  a  single  geometrical  parameter.  For 


simplicity,  we  assume  that  the  matrix  and  fiber  are  isotropic 


both  in  stiffness  and  thermal  expansion  coefficient.  Thus, 


‘=ljkf  °ljkt  “Ij 


IJki 


IJ  kt 


ik^Ji  it  kj 


(10) 


^IJki  “  ^  ^IJ^kt  ^^Ik^Jt  ^It^kJ^ 


(11) 


-  (a^  -  a)  Sj^jAT 


(12) 


In  the  above  equations,  6^^  is  the  Kronecker's  delta,  x  ( X^ ) 


and  y  (u  )  are  Lame's  constants  of  the  matrix  (fiber). 


a  (aj.)  is  the  thermal  expansion  coefficient  of  the  matrix 


(fiber)  and  aT  Is  the  change  in  temperature  (aT  >  0 


corresponds  to  the  temperature  drop). 


The  stress  field  in  the  fiber  is  obtained  from  Eqs.  (6), 


(8),  (9)  and  (10) 
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=  (1  -  Vf.){  X6j_.  +  2w 


»« 

®n>l 


(13) 


where  the  repeated  Index  Is  to  be  summed  over  1,  2,  and  3. 

«« 

After  solving  for  In  Eq.  (1)  by  use  of  Eqs.  (8)  and  (9), 
we  can  compute  the  stress  disturbance  o^j  In  the  fiber  from 
Eq.  (13). 


t  c 

B.  Yield  stresses  a  ,  o  and  work-hardening  rate 

y  y 


The  method  of  computing  the  yield  stresses  and  and 

y  y 

the  work-hardening  rate  of  the  composite  is  described 
briefly.  The  total  potential  energy  of  Fig.  3,  U,  is  given 


U  *  ?  /_  '^ij  '^ij  ■  ®ij  ■ 


a°jnj(u°  +  Uj^  +  u^)dS 


(14) 


o  ~ 


where  u^,  u^^  and  Uj^  are  the  displacement  components  corre¬ 
sponding  to  e°j »  and  e^^j,  respectively;  the  index  J 
preceded  by  a  comma  denotes  a  partial  differentiation  with 
respect  to  Xj ;  PI  is  the  boundary  of  D;  and  nj  is  the  J-th 
component  of  an  unit  vector  outer  normal  to  IDf.  Then,  the 

change  in  u  due  to  the  change  in  the  plastic  strain  6e?.  is 

J 

given  by 


^  V'- .  -•  “  * 


r. 


S.’ 


I  « 


)5efjdV 


(15) 


In  the  above  derivation,  the  Gauss’  divergence  theory  and  the 
following  equation  were  used 


(6U4  .  +  6u.  )dV  =  0 


'IJ 


1, 


Noting  that  the  plastic  strain  exists  only  In  the  matrix 
(D-G),  Eq.  (15)  Is  reduced  to 


6U 


^<3  -  Vijt 


(16) 


Under  the  uniaxial  stress  along  the  x^-axls  (o^),  and 


-  are  given  by 


'IJ 


IJ 


/^-l/2e. 


■IJ 


-l/2e. 


( 


(17) 


0 

0 

0 


where  the  six  components  of  and  e^.  are  expressed  In  the 

ij  ij 


order  of  (IJ)  *  11,  22,  33,  23,  31,  and  12,  and  Bp  Is  the 


plastic  strain  along  the  x^-axls.  On  the  other  hand,  the 
energy  dissipation  due  to  the  plastic  work  In  the  matrix, 
6Q,  Is  given  by 


fiQ  -  (l-V-)o„6e 


f '''y''''p 


(18) 
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where  o  is  the  flow  stress  of  the  matrix  for  the  bilinear 
model 


(19) 


where  and  Eip  are  the  Initial  yield  stress  and  tangent, 
modulus  of  the  matrix,  respectively,  and  e  Is  the  total 
strain.  Since  6U  +  6Q  =  0,  we  obtain 


(a33  -  0^^) 


(20) 


In  the  above  derivation,  =  22  “  ”  "I  ®p 

used.  Combining  the  solutions  of  In  Eq.  (13)  and  Eq. 
(20),  we  can  obtain  the  yield  stress  of  the  composite  In 
tension  (o^  =  a  )  and  that  In  compression  (a°  *  a  )  as 


''y'''’y  ■  '^0  *  ''i 


(21) 


■  Co  -  Cj  a^4T 


Similarly,  the  work-hardening  rate  of  the  composite,  Is 
obtained  as 


Et/Et  -  Cg  C3  (g^) 


(22) 


In  the  above  equations,  Cf  ,  Cj^,  C2,  and  are  functions  of 


the  mechanical  properties  of  the  matrix  and  fiber  and  the 
fiber  aspect  ratio  £/d,  E  is  the  Young's  modulus  of  the 
matrix. 

C.  Thermal  Residual  Stress 


When  the  composite  Is  cooled  down  by  aT,  the  thermal 
residual  stress  Is  Induced  In  the  composite.  The  theoretical 
model  for  this  problem  Is  the  same  as  that  shown  In  Pig.  ^ 
except  that  =  e^^  =  0  In  the  present  case.  Hence,  the 
formulation  up  to  Eq.  (13)  Is  valid  and  will  be  used  to 
compute  the  thermal  residual  stress.  Once  the  stress  within 
the  fiber  Is  computed,  the  tresses  Just  outside  the 

fiber  can  be  obtained  by  the  following  relatlon^^^^ 


(out)  ^  „(ln)  „  r  p  »» 

pq  pq  ^pqmn^  ^k£lj®lj”£% 


(X-^2u)6j^m~ 

VI  ( X  +  2y  ) 


(23) 


where  Is  given  by  Eq.  (13),  e??  Is  solved  by  Eq.  (1), 

pq 

X  and  g  are  previously  defined  and  n^^  Is  the  1-th  component 
of  the  unit  vector  prependicular  to  the  surface  of  the  fiber. 
If  one  wants  to  compute  the  stresses  Just  outside  the  equator 
of  the  ellipsoidal  fiber  where  n  “  (1,  0,  0),  we  can  obtain 
the  stresses  In  polar  coordinates,  o  ,  o.  and  o  there  as 

P  w  Z 


(out)  ^  (In) 


(out) 


^  2-It^  eji  . 


V  ^33' 


(2H) 


(out)  ^  (in)  ^  2y{^  e**  +  ^  e**} 
z  z  ^‘1-v  11  1-v  33' 


where  v  Is  the  Poisson’s  ratio  of  the  matrix. 

Next  the  average  stress  field  In  the  matrix, 
be  also  obtained  from  Eqs.  (3),  (8)  and  (9)  as 


<a,  .>  = 

ij  ra 


^f  ^IJkt ^^kimn^mn  ®ki ^ 


(25) 


The  average  stress  In  the  matrix,  ^°z^m* 

will  be  computed  by  setting  iJ  *  11,  22,  and  33  In  Eq.  (25). 
It  should  be  noted  here  that  <o>^  = 

assumptions  of  an  aligned  short  fiber  composite  which  results 
In  the  transverse  Isotropy  of  the  volume  average  quantity. 
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5.  Theoretical  Results 

We  vflll  apply  the  present  theoretical  model  to  two  types 
of  MMCs,  short  fiber  MMC  (SIC  whlsker/606l  Al)  and  continuous 
fiber  MMC  (AI2O3-Y/5056  Al). 

A.  Short  Fiber  MMC 

The  thermo-mechanical  data  of  the  matrix  and  whisker  for 
the  theoretical  calculations  are  obtained  from  the  stress- 
strain  curve  of  the  matrix  (Pig.  2a)  and  the  material 
properties  handbook. 

Annealed  6061  Al  matrix  (bilinear  model): 

E  =  l<7.5  GPa 
Eip  *  2.3  GPa 

a°  =  JI7.5  MPa  (26) 

y 

V  =  0.33 

a  =  23.6  X  lO'^/K 

SIC  Whisker: 

E^  »  427  GPa 

V  -  0.17  (27) 

w 

«  4.3  X  lO'^/K 

£/d  -  1.8 


where  the  average  value  of  the  fiber  aspect  ratio  (l/d)  was 
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and  the  bilinear  stress  strain  curve  of  the  matrix 
Is  Indicated  by  a  dash-dot  line  in  Fig.  2a.  The  temperature 
drop  AT  Is  defined  as 


AT  =  Tj  -  Tq 


(28) 


where  T^  Is  taken  as  the  temperature  below  wh  ch  dislocation 
generation  Is  minimal  during  the  cooling  process and  Tq  Is 
the  room  temperature.  Thus,  for  the  present  composite  system 
at  Is  set  equal  to  200  K. 

c  t 

The  theoretically  predicted  o  and  a  can  be  obtained 

y  y 

from  Eq.  (21).  The  differences  between  and  (a^  -  o^) 

y  y  y  y 

for  various  Vj.  are  shown  as  the  solid  line  in  Fig.  5*  The 

c  t 

differences  between  a  and  o  Increase  with  Increasing  Vj.. 

y  y 

c  t 

The  experimentally  determined  differences  between  a  and  a 
are  represented  by  the  open  circles  In  Pig.  5.  The  experi¬ 
mental  values  of  the  yield  stresses  are  taken  as  the  stresses 
0.2%  off-set  strain.  It  follows  (Pig.  5)  that  good  agreement 
Is  obtained  between  the  experimental  and  theoretical  results 
of  the  difference  In  the  compressive  tensile  yield  stresses. 
However,  the  experimental  values  of  are  greater  than  the 

y 

theoretical  values  of  a°,  and  similarly  the  experimental 
values  of  are  greater  than  the  theoretical  values  of  o^. 

Prom  the  data  given  by  Eqs.  (26)  -  (27)  and  the  use  of 
Eq.  (25),  we  have  computed  the  stresses,  and  Its  average 
value  In  the  matrix  and  plotted  schematically  the  a.  (solid 

ii 

curve)  and  <a- >  (dashed  curve)  along  the  x,  and  x,  axes  In 
li  m  ^  j 


■li 
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I 


-v« 

■'v 
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Pig.  6. 

Next,  the  thermal  residual  stresses  averaged  In  the 
matrix  of  the  SIC  whlsker/606l  A1  are  predicted  by  Eq.  (25) 
and  the  results  on  <o„>„  and  <0.  >  are  plo1|fced  In  Fig.  7  as  a 
function  of  the  volume  fraction  of  whisker  (Vf),  where  the 
subscripts,  T  and  L  denote  the  component  along  the  transverse 
direction  (r  and  e)  and  longitudinal  direction  (z).  The 
average  theoretical  thermal-residual  stress  Is  predicted  to 
be  tensile  In  nature.  This  theoretical  prediction  Is  In 
general  agreement  with  the  experimental  results  (Table  I). 
Also,  the  residual  stress  is  experimentally  observed  to  be 
larger  In  the  transverse  direction  and  a  larger  transverse 
stress  Is  shown  In  Fig.  7. 

For  the  case  of  spherical  SIC  composites  the  procedure 
for  obtaining  the  difference  In  0°  and  Is  Identical  to  the 

y  y 

whisker  case,  with  one  exception,  l.e.,  £/d  =  1  In  the 

spherical  case.  The  theoretically  predicted  difference 
c  t 

between  a  and  o  for  Vg  =  0,2  Is  zero,  but  experimentally  it 

y  y 

c  t 

was  observed  that  a  was  slightly  greater  than  a  by  ~  13 

MPa. 

B.  Continuous  Fiber  MMC 

The  present  model  Is  also  applied  to  continuous  fiber 
(AI2O2-Y)  reinforced  5056  A1  composite  with  Vj.  ■  0.5  (Al202“ 
y/5056  A1  In  short  abbreviation)  to  predict  the  yield  stress 
In  tension  and  compression,  and  also  the  work-hardening  rate. 
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The  material  properties  of  the  matrix  metal  and  fiber 

are: 

Annealed  5056  A1  matrix  (bilinear  model): 

E  =  72. GPa 
E>p  =  6.25  GPa 


a°  =  152  MPa 
y 

(29) 

V  =  0.33 

a  =  23.6  X  10“^/K 

AI2O3-Y  Flber^^^^ 

Ef.  =  210  GPa 

Vj.  =  0.2 

(30) 

0-  »  8.8  X  10"^/K 

In  this  computation  AT  =  200  K  was  used. 

Figure  8  shows  a  comparison  between  the  experimental 
results^^^^  (solid  curve  and  filled  circle)  and  the 
theoretical  ones  (dashed  curve  and  open  circle)  where  both 
filled  and  open  circles  denote  the  yield  stresses.  It 
follows  from  Fig.  8  that  the  yield  stress  and  the  work- 
hardening  rate  of  the  compressive  stress-strain  curve  are 
well  predicted  by  the  present  model,  while  for  tho  tensile 
yield  stress,  the  present  model  underestimates  the 
experimental  results. 


The  theoretical  predictions  and  experimental  results  are 
in  very  good  agreement  in  most  cases. 

The  theoretical  model  is  based  on  an  extension  of  pre¬ 
vious  work  by  Eshelby^^^\  Mura  and  Taya^^^^,  Tanaka  and 
Mori and  Wkashlma,  et  al.^^®^  and  it  can  predict  the 
yield  stress  in  tension  and  compression  and  the  thermal 
residual  stress.  The  predicted  values  of  the  yield  stress  of 
the  whisker  composite  in  tension  and  compression  are  less 
than  the  experimentally  determined  values  of  the  yield  stress 
In  tension  and  compression.  This  difference  between  the 
theoretical  and  experimental  stresses  Is  due  to  the  fact  that 
the  *  0  curves  of  stress  vs.  strain  were  used  for  the 
matrix  In  the  composite  cases.  It  has  been  demonstrated  by 
Arsenault  and  Fisher ^®^  and  Vogelsang,  et  al.^^^  that  there 
is  a  much  higher  dislocation  density  In  the  annealed 
composite  matrix  than  In  the  «  0  material.  Therefore,  the 
matrix  Is  stronger  than  the  annealed  Vj.  »  0  matrix  alloy. 
However,  this  Increase  In  matrix  strength  due  to  a  higher 
dislocation  density  does  not  influence  the  difference  In 
yield  stress  In  the  compression  vs.  tension.  The  theoretical 
prediction  Is  that  the  yield  stress  of  the  whisker  composite 
In  compression  should  be  higher  than  the  yield  stress  In 
tension,  and  this  Is  exactly  what  Is  observed  experimentally, 
due  to  an  average  residual  tension  stress  In  the  sample. 

This  residual  stress  reduces  the  yield  stress  In  tension  and 
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increases  the  yield  stress  In  compression. 

The  theoretical  prediction  that  the  compressive  yield 
stress  Is  higher  than  the  tensile  yield  stress  suggests  that 
a  tensile  residual  stress  exists  In  the  matrix  of  composite. 
An  average  tensile  residual  stress  is  predicted  In  the  matrix 
as  shown  In  Figs.  6  and  7. 

Though  the  present  model  Is  based  on  the  assumption  that 
all  fibers  are  aligned  along  a  specified  direction  (the  x^- 
axis  In  Pig.  3),  SIC  whiskers  In  SIC  whlsker/6061  A1  compos¬ 
ites  are  mlsorlented.  The  mlsorlentatlon  of  these  whiskers 
Is  observed  to  be  enhanced  as  Vf  decreases.  Thus  the  values 
of  and  a®  predicted  by  the  present  model  tends  to  over- 

y  y 

estimate  the  actual  values.  The  fiber  aspect  ratio  (l/d)  of 
SIC  whiskers  has  been  observed  to  be  variable l.e., 

0.5  <  l/d  <  25  with  most  of  A/d  being  around  1.  In  the 
present  model  we  have  used  the  average  value  of  £/d,  1.8  to 
predict  the  yield  stress  and  work-hardening  rate  of  the 
composite.  Takao  and  Taya^^^^  have  recently  examined  the 
effect  of  variable  fiber  aspect  ratio  on  the  stiffness  and 
thermal  expansion  coefficients  of  a  short  fiber  composite  and 
concluded  that  the  results  predicted  by  using  the  average 
value  of  A/d  are  not  much  different  from  those  of  using 
actual  variable  aspect  ratio  unless  the  distribution  of  A/d 
Is  extraordinaryly  scattered.  Thus,  the  use  of  the  average 
value  of  A/d  In  the  present  model  can  be  Justified. 


22 


The  present  model  Is  strongly  dependent  on  AT  (■  T^  - 
Tq)  as  it  appears  in  Eqs.  12,  21  and  22.  A  choice  of  Tj^  was 
made  by  noting  the  movement  of  dislocation  generation  being 
minimal^ as  discussed  before.  If  temperature  Tj^  can  be 
continued  to  correspond  to  the  temperature  Tp  above  which  the 
relaxation  of  Internal  stress  Is  highly  promoted  and  below 
which  the  deformation  of  the  matrix  metal  Is  basically 
controlled  by  elastic  and  time-independent  plastic  behavior, 
then  the  present  model  which  Is  based  on  an  elastlc/time- 
Independent  plasticity  becomes  valid.  If  T2  la  not 
approximately  equal  to  Tp,  a  more  refined  model  to  account 
for  elastlc/tlme-lndependent  plastic/creep  behavior  of  the 
matrix  metal  during  a  temperature  excursion  from  T^  to  Tq 
must  be  employed.  Such  a  model  for  one-dlraensional  stress 
analysis  has  been  proposed  by  Garmong^^^^. 


7.  Conclusions 


It  is  possible  to  arrive  at  the  following  conclusions 
from  the  theoretical  model  of  the  thermal  residual  stresses 
and  the  experimental  results. 

1.  The  theoretical  model  predicts  a  higher  compression  yield 
stress  than  tensile  yield  stress  for  the  whisker 
composite,  which  Is  In  agreement  with  the  experimental 
results.  The  higher  compressive  yield  stress  Is  due  to  a 
tensile  residual  stress.  This  residual  stress  Is  due  to 
the  difference  In  coefficient  of  thermal  expansion 
between  fiber  and  matrix. 

2.  The  absolute  magnitude  of  the  predicted  yield  stress 
(tensile  and  compressive)  Is  less  than  the  experimental 
yield  stress  due  to  an  Increase  In  matrix  strengthening. 
The  Increased  matrix  strength  Is  due  to  a  higher 
dislocation  density  In  the  matrix;  the  higher  dislocation 
density  Is  a  consequence  of  the  difference  In  coefficient 
of  thermal  expansion  between  the  A1  alloy  matrix  and  the 
SIC. 

3.  The  residual  stress.  In  the  region  between  SIC  whiskers. 
Is  tensile  and  this  Is  the  likely  region  where  plastic 
deformation  would  begin.  For  In  this  region  the  matrix 
contains  a  lower  dislocation  density  than  adjacent  to  the 
SIC  whisker,  l.e..  In  this  region  the  matrix  Is  weaker. 


The  theoretical  model  predicts  that  the  compressive  yield 
stress  should  equal  the  tensile  yield  stress  for 
composites  containing  spherical  SIC.  Experimentally  It 
Is  observed  that  the  tensile  yield  stress  Is  slightly 
larger  than  the  compressive  yield  stress. 

5.  The  present  model  can  predict  the  yelld  stress  and  work¬ 
hardening  rate  of  continuous  flber/5056  A1 

composite  reasonably  well  except  for  o^. 

6.  In  the  short  fiber  case,  the  theoretical  model  predicts 
that  the  matrix  should  be  In  tension  and  experimental  X- 
ray  data  confirms  this  prediction. 
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Table  I 


Thermal  Residual  Stress  (Tensile) 
X-Ray  Measurements 


Material 

Transverse 

Longitudinal 

MPA 

MPA 

0  V%  Whisker  SIC 

6061  Matrix 

0.0 

0.0 

5  V*  Whisker  SIC 

6061  Matrix 

i»08 

35 

20  V%  Whisker  SIC 

6061  Matrix 

231 

58 

Wrought 

1100  A1 

0.0 

0.0 
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Figure  Captions 

Figure  1  A  schematic  diagram  of  the  sample  configuration 
which  was  used  In  the  tension  and  compression 
testing.  The  dimensions  are  In  millimeters. 

Figure  2a  The  absolute  stress  vs.  strain  curves  for  tension 
and  compression  test  of  0  volume  %  whisker 
material  In  the  annealed  condition. 

Figure  2b  The  absolute  value  of  stress  vs.  strain  for  5 
volume  %  whisker  composite  In  the  annealed 
condition. 

Figure  2c  The  absolute  value  of  the  stress  as  a  function  of 
strain  for  20  volume  %  whisker  material  In  the 
annealed  condition. 

Figure  3  Theoretical  model;  actual  case. 

Figure  4  The  equivalent  Inclusion  model  converted  from 
Figure  3. 

Figure  5  The  difference  In  yield  stress  between  compression 
and  tensln  as  a  function  of  volume  %  silicon 
carbide  whisker. 

Figure  6  A  schematic  of  the  stress  distribution  In  the 
matrix  and  In  the  reinforcement  due  to  the 
difference  In  thermal  coefficient  of  expansion 
between  the  slllcon-carblde  and  aluminum. 


Figure  7  Predicted  residual  stress  In  the  matrix  for  the 
transverse  and  longitudinal  directions. 
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Figure  8  Tension  and  compression  stress-strain  curves  of 

continuous  AI2O2-Y  flber/5056  A1  composite  with  Vj. 
=  The  experimental  and  theoretical 

results  are  denoted  by  solid  dashed  curves, 
respectively.  The  filled  OAd  open  circles  denote 
the  yield  stress  of  the  experimental  and 
theoretical  results,  respectively. 
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which  was  used  in  the  tension  and  compression 
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Figure  2a  The  absolute  stress  vs.  strain  curves  for  tension 


and  compression  test  of  0  volume  %  whisker 
material  In  the  annealed  condition. 


Figure  2b  The  absolute  value  of  stress  vs,  strain  for  5 
volume  %  whisker  composite  In  the  annealed 


condition 


Figure  2c  The  absolute  value  of  the  stress  as  a  function  of 
strain  for  20  volume  %  whisker  material  in  the 
annealed  condition. 
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Figure  6 


A  schematic  of  the  stress  distribution  In  the 
matrix  and  In  the  reinforcement  due  to  the 
difference  In  thermal  coefficient  of  expansion 
between  the  slllcon-carblde  and  aluminum. 
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Figure  8 


Tension  and  compression  stress-strain  curves  of 
continuous  Al202"Y  fiber/5056  A1  composite  with 
-  0.5(26).  The  experimental  and  theoretical 
results  are  denoted  by  solid  dashed  curves, 
respectively.  The  filled  9ad  open  circles  denote 
the  yield  stress  of  the  experimental  and 
theoretical  results,  respectively. 


